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1)  Abstract 


Smart  composites  comprising  of  electro-active  components  are  appealing  for  multifunctional 
systems  with  the  ability  to  respond  to  various  external  stimuli,  such  as  morphing  structures. 
Piezocomposites  having  piezoelectric  ceramics  and  polymers  are  promising  for  morphing 
structures.  The  electro-mechanical  performances  of  piezoelectric  ceramics  and  polymers  are 
nonlinear  and  time-dependent,  which  are  significantly  affected  by  temperatures.  The  differences 
in  the  responses  of  piezoelectric  ceramics  and  polymers,  together  with  various  microstructural 
morphologies  will  have  significant  effects  on  the  overall  performance  of  multifunctional 
morphing  systems. 

This  project  is  aimed  at  understanding  the  effect  of  coupled  thermal,  electrical,  and  mechanical 
responses,  including  loading  rate  (time)  effect,  of  the  constituents  (piezoelectric  ceramics  and 
polymers)  on  the  overall  multi-field  responses  of  electro-active  composites.  For  this  purpose,  a 
multi-scale  framework  that  consists  of  constitutive  models  for  the  constituents  incorporating 
field  coupling  and  time-dependent  effects  and  heat  conduction;  nonlinear  micromechanical 
models  for  various  composite  reinforcements  including  functionally  graded  composites;  and 
large  scale  multifunctional  structural  analyses  is  formulated.  Limited  experiments  on 
piezocomposites  are  conducted  in  order  to  provide  useful  information  on  the  effect  of 
microstructural  morphologies  on  the  overall  multi-field  response  of  composites.  The  constitutive 
and  micromechanical  models  are  capable  in  capturing  the  time-dependent  multi-field  responses 
of  the  smart  composites.  The  multi-scale  framework  can  predict  the  overall  performance  and 
perform  active  controls  in  order  to  achieve  desired  performance  of  the  structures  under  various 
stimuli.  For  precisely  controlling  the  shapes  of  structures,  it  is  necessary  to  consider  the  different 
time-dependent  behaviors  of  the  constituents  and  microstructural  morphologies  in  the  active 
composites. 
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3)  Statement  of  Problem  Studied 


This  project  is  aimed  at  understanding  the  effect  of  coupled  thermal,  electrical,  and  mechanical 
responses,  including  loading  rate  effect,  of  the  constituents  on  the  overall  multi-field  responses 
of  electro-active  composites.  Laminated  composites,  having  integrated  piezocomposite 
transducers  of  piezoelectric  ceramics  fiber  or  particle  and  homogeneous  matrix  (polymers  or 
metals),  suitable  for  multifunctional  morphing  structures  are  studied.  The  ultimate  goal  is  to 
formulate  a  multi-scale  framework  that  consists  of  constitutive  models  for  constituents 
incorporating  coupled  deformation,  heat  conduction,  and  piezoelectric  analyses  with  stress, 
temperature,  electric  field,  and  time  dependent  responses;  nonlinear  micromechanical  models  for 
various  composite  reinforcements  including  functionally  graded  composites;  and  large  scale 
multifunctional  structural  analyses.  In  addition  to  the  multi-scale  development,  limited 
experiments  on  piezocomposites  are  conducted  in  order  to  study  micro-scale  response  of  the 
composite  and  provide  useful  insight  into  the  effect  of  microstructural  morphologies  on  the 
overall  (macroscopic)  multi-field  response  of  composites. 

The  project  consists  of  the  following  components: 

1)  Development  of  nonlinear  time-dependent  electro-mechanical  constitutive  models  for  the 
active  constituent,  i.e.,  piezoelectric  ceramics,  and  nonlinear  viscoelastic  and/or 
viscoplastic  models  that  incorporate  the  effect  of  temperatures  and  energy  dissipations, 
which  is  used  for  the  polymeric  constituent.  The  constitutive  models  are  derived  based  on 
a  phenomenological  approach  within  continuum  mechanics. 

2)  Formulation  of  micromechanics  models  for  fiber  and  particle  composites.  Different 
microstructural  morphologies  of  composites  are  considered,  which  are  simplified  unit¬ 
cell  models,  composites  with  uniformly  and  randomly  distributed  inclusions  (long  fibers 
and  solid  spheres),  and  composites  with  microstructural  details  obtained  from  the 
scanning  electron  microscope.  The  micromechanical  models  are  used  to  study  the 
nonlinear  coupled  heat  conduction  and  deformation  in  viscoelastic  composites  and 
nonlinear  time-dependent  thermo-electro-mechanical  response  of  active  composites. 

3)  Implementation  of  the  constitutive  models  and  micromechanical  models  within  finite 
element  (FE)  and  finite  difference  (FD).  FE  method  is  then  used  for  analyzing  smart 
composite  structures/structural  components  undergoing  coupled  thermo-electro- 
mechanical  stimuli,  while  incorporating  detailed  nonlinear  and  time-dependent  responses 
of  the  constituents  within  the  smart  composite  structures. 


4 


4)  Summary  of  the  Most  Important  Results 


This  project  presents  multi-scale  analyses  of  smart  composite  structures  that  consist  of  electro¬ 
active  materials  and  non-active  components.  The  focus  is  on  studying  flexible  smart  composites 
with  the  ability  to  undergo  shape  changes  in  responses  to  external  stimuli.  The  electro-active 
materials  considered  here  are  piezoelectric  ceramics  and  piezoelectric  composites  that  comprise 
of  piezoelectric  ceramics  and  ductile  matrix,  such  as  polymers  and  metals.  The  non-active 
structural  components  are  made  of  fiber  or  particle  reinforced  composites.  When  subjected  to 
mechanical  and  non-mechanical  (thermal  and  electrical)  stimuli,  the  smart  composite  structures 
would  experience  nonlinear  and  time-dependent  coupling  responses  due  to  the  different 
behaviors  of  the  constituents.  This  project  aims  at  integrating  different  nonlinear  and  time- 
dependent  responses  of  the  constituents  in  smart  composites  in  order  to  predict  the  overall 
performance  of  smart  structures. 


4.1)  Experimental  Observation  on  the  Nonlinear  Time-dependent  Response  of 
Piezocomposites 

Piezocomposites,  comprising  of  long  unidirectional  piezoelectric  ceramics  fibers  embedded  in 
polymeric  matrix,  also  known  as  active  fiber  composites  (AFCs)  are  studied  in  terms  of  their 
electro-mechanical  behaviors  at  different  temperatures.  The  studied  AFC  specimens,  which  were 
purchased  from  Advanced  Cerametrics  Incorporated,  consist  of  PZT  5  A  fibers  dispersed  in  an 
epoxy  matrix.  Kapton  layers  and  electrode  fingers  are  placed  on  both  surfaces  of  the  specimens, 
resulting  in  electric  field  along  the  longitudinal  fiber  axis.  The  AFC’s  concept  was  developed  by 
Bent  (1994)  and  Hagood  et  al.  (1993),  which  is  illustrated  in  Fig.  1  (left).  A  photograph  of  the 
tested  specimen  is  shown  in  Figure  1  (right). 


Figure  1  Microstructural  representation  of  AFC  (left)  and  fracture  of  an  AFC  specimen  after  tensile  test 
(right) 


The  experimental  tests  are  conducted  in  order  to  understand  the  quasi-static  and  hysteretic 
mechanical  and  electromechanical  responses  of  AFC  at  various  temperatures  and  loading  rates 
(see  Ben-Attitalah  et  al.  2014  for  a  discussion  on  experimental  tests  and  results).  First,  the 
mechanical  properties  of  AFCs  are  evaluated  by  applying  tensile  tests  in  the  longitudinal  fiber 
direction,  at  different  rates  and  temperatures:  25°C,  50°C,  and  75°C.  Figure  2  shows  examples  of 
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the  uniaxial  stress-strain  responses  and  Table  1  present  a  summary  of  mechanical  responses, 
which  indicate  rate  and  temperature  dependent  behaviors  of  AFCs.  The  specimens  are  also  tested 
under  creep  loadings  at  different  stresses  and  temperatures,  and  under  relaxations  at  various 
strains  and  temperatures.  The  results  show  significant  time-dependent  behaviors.  Hysteretic 
responses  of  AFCs  are  also  studied  under  various  amplitude  and  frequencies  of  electric  field 
input  and  at  different  temperatures.  Figure  3  illustrates  hysteretic  polarization  response  of  AFCs 
at  different  frequencies  and  temperatures.  The  amplitude  of  electric  field  is  taken  as  25  kV/cm. 
The  responses  are  compared  to  those  of  pure  PZT5A  specimens,  which  shown  quite  different 
behaviors.  Figure  4  presents  the  responses  of  AFCs  at  higher  amplitude  of  electric  field,  which  is 
40  kV/cm.  In  Figs.  3  and  4,  high  electric  fields  are  applied  so  that  the  piezoelectric  ceramics 
experience  polarization  switching.  When  relatively  low  electric  field  inputs  are  prescribed  on 
polarized  AFCs,  hysteretic  electro-mechanical  responses  are  observed,  as  shown  in  Fig.  5.  The 
tests  are  also  done  at  different  temperatures  and  frequencies.  It  is  seen  that  the  mechanical, 
electrical,  and  electro-mechanical  responses  of  AFCs  depend  strongly  on  rate  of  loading, 
temperatures,  and  magnitude  of  the  mechanical  and  electric  field  inputs.  Thus,  in  order  to  better 
understand  the  overall  performance  of  active  composite  structures  made  of  electro-active 
materials,  it  is  necessary  to  include  the  knowledge  of  time-dependent  and  field  coupling  effects 
of  the  different  constituents  in  the  active  structures. 


Figure  2  Uniaxial  stress-strain  under  different  rates  at  room  temperature  (left)  and  at  different  temperature 
and  loading  rate  25min/min  (right) 


Table  1  Mechanical  response  of  piezocomposites  at  various  temperatures  and  high  strain  rate 
(average  from  3  repeated  tests) 


Displ.  Rate 
(mm/min) 

Elastic  modulus  (GPa) 
25°C  50°C  75°C 

Tensile  strength  (MPa) 
25°C  50°C  75°C 

Max  strain  (%) 

25°C  50°C  75°C 

0.5 

9.78 

6.25 

5.96 

51.29 

44.04 

39.90 

17.94 

30.94 

34.13 

5.0 

16.58 

7.66 

4.05 

43.42 

51.13 

44.50 

5.37 

21.60 

26.24 

25.0 

16.47 

8.49 

6.27 

44.68 

47.03 

39.15 

4.93 

23.80 

36.78 
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Figure  3  Hysteresis  response  of  PZT  and  active  fiber  composites  at  various  frequencies  at  temperatures 
25‘>C  and  75“C. 


Ferroelectric  hysteresis  loop  for  an  AFC  4kVpp 


Figure  4  Hysteretic  polarization  at  room  temperature  with  maximum  eleetrie  field  40  kV/em 
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Figure  5  Electro-mechanical  responses  of  polarized  AFCs  at  frequency  IHz  and  different  amplitude  of 

electric  fields. 


4.2)  Nonlinear  Time-dependent  Constitutive  Models 

Polymers  and  piezoelectric  ceramics  of  perovskite  structures  are  considered  in  this  project,  and 
they  are  shown  to  experience  nonlinear  and  time-dependent  responses  under  mechanical  and 
non-mechanical  inputs,  as  discussed  in  the  experimental  tests  above  and  in  several  literature. 
Nonlinear  time-dependent  constitutive  models  incorporating  field  coupling  effects  are  developed 
for  different  constituents  in  the  active  composite  structures.  The  multi-field  responses  of  the 
constituents  are  assumed  to  follow  general  integral  functions  in  order  to  incorporate  histories  of 
loading. 


4.2.1  Thermo-mechanical  models  of  polymers 


In  case  of  polymers,  a  nonlinear  single  integral  model  of  Schapery  (1969)  is  used  to  represent  the 
viscoelastic  response  of  polymers.  The  model  also  include  the  effect  of  temperatures  and 
dissipation  of  energy  that  is  converted  into  heat.  The  model  is  then  extended  in  order  to  include 
the  viscoplastic  response  of  polymers.  Linearized  strain  measures  are  considered  in  the  model. 


which  are  e,,  =  — 
2 


(w.  .  +w . .) 


,  where  .  is  the  displacement  gradient.  The  nonlinear  thermo¬ 


viscoelastic  model  is  given  as  (Muliana  and  Khan  2008): 


4  +^i 


(1) 


where  the  superscript  t  denotes  the  corresponding  variables  at  time  t,  Do  and  dD  are  the 
instantaneous  elastic  and  transient  compliances  under  a  uniaxial  (extensional)  creep  loading, 
respectively.  The  nonlinearity  is  incorporated  through  parameters  go,  gi,  g2  which  are  functions 
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of  current  temperature  T  ‘  and  the  effective  stress  cr  ‘ .  The  reduced  time  xf/‘  is  defined  as 
y/*  =\i/U^=  [  where  the  functions  ar  and  are  the  temperature-  and  stress  dependent 

0  ^T^cj 

parameters.  Using  the  energy  equation  that  incorporates  the  viscoelastic  dissipation,  and  stress- 
and  temperature  dependent  material  parameters,  thermo-mechanical  coupling  becomes: 


pc^{r)t  +  T‘a(T‘)&‘^+r 


Sgo(o-,T‘)  ^(0)-,, 


dT 


-Tc7i,+r 


dT 


-r&^+r 


dT 


ijkl  ij  kl 


0  rn=\ 


^  d  ‘ 


J dDyti  (¥' 


where  a,y(r^)  are  the  components  of  coefficient  of  thermal  expansion  (CTE)  tensor  and  can 

depend  on  temperatures,  /7is  the  density,  c^{T‘)  is  the  temperature  dependent  specific  heat  at  a 

constant  stress,  and  qt  is  the  scalar  component  of  the  heat  flux.  The  above  models  are 
implemented  within  finite  element,  which  are  used  to  analyze  coupled  thermo-viscoelastic 
response  of  polymeric  structures.  The  models  are  also  integrated  to  micromechanics  models  and 
used  to  study  fully  coupled  heat  conduction  and  thermo-mechanical  responses  of  composite 
structures.  Detailed  discussion  on  the  model  and  numerical  implementation  are  given  in  Khan 
and  Muliana  (2012a  and  b). 

Figure  6  illustrates  comparisons  between  experimental  data  and  finite  element  analyses  of  a 
torsional  polymeric  bar  under  oscillations  at  different  frequencies.  The  adiabatic  heating  is 
considered.  It  is  seen  that  the  temperature-dependent  material  parameters  affect  the  heat 
generation  in  the  polymeric  bar  due  to  the  viscoelastic  dissipation.  Another  example  is  shown  on 
a  simply  supported  beam  subjected  to  a  uniformly  distributed  load  with  a  sinusoidal  function. 
Three  different  materials  are  considered  for  the  beam,  which  are  pure  polymers,  particle 
reinforced  polymers  with  a  uniform  distribution  of  particulate  inclusions,  and  functionally  graded 
materials  with  gradual  variations  of  the  particles  in  polymeric  matrix.  Figure  7  depicts 
temperature  fields  in  the  beam  at  two  different  times.  Highest  temperatures  are  generated  in  the 
polymeric  beam,  as  expected,  since  the  unreinforced  polymer  experiences  a  more  pronounced 
viscoelastic  behavior.  Adding  elastic  particles  to  the  polymers  reduces  the  overall  viscoelastic 
responses.  Gradually  varying  the  particle  compositions  can  minimize  the  heat  generations. 
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Experiment  (Tauchert,  1967a) 

Proposed  Model  (Adiabatic  Condition) 
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0=  Bq  sinoot 
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3: 
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^  ^ _ _ 


1000  1500 
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Figure  6  Comparison  of  experimental  and  finite  element  results  under  torsional  oscillations  at  different 
constant  frequencies  and  angle  of  twist  of  6o=30°:  a)  Adiabatic  solution  with  temperature  independent 
material  properties,  b)  Adiabatic  solution  with  temperature-dependent  material  properties. 


Time=12.5  mins 


Time=25  mins 


304.5  K  ■  306  K 

Figure  7  Temperature  fields  due  to  the  dissipation  effects  during  cyclic  bending  in  viscoelastic  beams 
comprising  of  a)  homogeneous,  b)  composite  (12.5%  particle  contents),  and  c)  functionally  graded 
material 


303  K 


Some  polymers  experience  inelastic  deformations  when  subjected  to  relatively  large  mechanical 
loadings,  resulting  in  permanent  deformations.  In  this  study,  the  effect  of  inelastic  deformations 
on  the  time-dependent  behavior  is  incorporated  by  including  the  viscoplastic  component.  The 
total  strain  is  additively  decomposed  into  the  viscoelastic  and  viscoplastic  parts.  The  single 
integral  model  in  Eq.  (1)  is  used  for  the  viscoelastic  part,  while  the  viscoplastic  part  follows  the 
Perzyna  (1971)  model: 

S,J  (3) 

where  sjf’’  is  the  viscoplastic  strain  rate  at  current  time  t  and  V  is  the  magnitude  of  the 
viscoplastic  strain  rate.  The  stress  dependent  yield  function  F(a^/c')  is  expressed  in  terms  of  the 
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effective  stress 


<7‘  = 


S‘.S‘.  , 

V  V  ’ 


where  5'.  is  the  deviatoric  stress  components  and  the  accumulated 


viscoplastic  strain  is  /r'  =  |  k'ds',  id  = 

0 

the  yield  surface  F  is  derived  as  dF  /  da.j  which  is  the  direction  of  the  viscoplastic  strain  rate. 

Another  viscoplastic  model  that  has  been  considered  is  an  endochronic  viscoplastic  model  of 
Valanis  (1971),  which  is  used  when  the  inelastic  deformation  occurs  at  relatively  low  stresses. 
The  constitutive  models  are  implemented  in  finite  elements,  which  are  discussed  in  Kim  and 
Muliana  (2010).  Figures  8  and  9  illustrate  the  numerical  results  of  the  viscoelastic-viscoplastic 
models  for  polyethylene.  Experimental  data  are  obtained  from  Lai  and  Bakker  (1995).  The 
viscoelastic-viscoplastic  model  is  integrated  to  the  micromechanics  models  of  fiber  and  particle 
reinforced  polymers,  which  is  discussed  in  Section  4.3. 


.  The  normal  direction  at  the  stress  point  on 


ij 


Figure  8  Strain  responses  under  two  strep  loadings 


Strain  (mm/nim) 

Figure  9  Stress-strain  under  different  rates 
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4.2.2  Electro-mechanical  models  for  piezoelectric  ceramics 

The  piezoelectric  ceramics  considered  here  are  of  perovskite  structures.  Each  crystallite  within 
the  polycrystalline  structures  can  undergo  spontaneous  polarization,  which  is  categorized  as 
ferroelectric  materials  (see  Lines  and  Glass  2009).  Several  available  experimental  data  on  lead 
zirconate  titanate  (PZT)  and  barium  titante  shows  nonlinear  and  time-dependent  electro¬ 
mechanical  responses,  e.g.,  Crawley  and  Anderson  (1990),  Fang  and  Li  (1999),  Zhou  and 
Kamlah  (2006).  The  nonlinearity  becomes  significant  under  high  electric  field  inputs,  and  the 
piezoelectric  ceramics  experience  polarization  switching  when  high  electric  fields  above  the 
coercive  limit  are  prescribed.  Nonlinear  electro-mechanical  constitutive  models  for  piezoelectric 
ceramics  incorporating  time-(rate-)dependent  effect  are  formulated  for  simulating  responses  of 
piezoelectric  ceramics  under  various  histories  of  electro-mechanical  stimuli. 


The  rate-dependent  polarization  switching  constitutive  model  for  ferroelectric  ceramics 
presented  by  Sohrabi  and  Muliana  (2013)  is  written  as: 

(0  =  ■SjX/  +  .  (■») 

(5) 

where  the  superscript  t  denotes  the  current  time  variable.  s\. ,  cr'. ,  D\  and  P'  are  the 

components  of  the  strain,  stress,  electric  displacement  and  polarization  tensors  at  current  time  t, 
respectively.  The  material  properties  are  the  elastic  compliance  tensor  determined  at  a 

constant  electric  field  and  dielectric  coefficient  tensor  k.j  calibrated  at  a  constant  stress.  Here 
is  the  scalar  component  of  the  piezoelectric  coefficient,  which  depends  on  current  polarization 
Pj ;  X3  direction  is  chosen  as  the  poling  axis.  Different  functions  can  be  chosen  for 

depending  on  the  experimental  data.  The  polarization  switching  is  due  to  an  application  of  cyclic 
electric  field  £[ .  The  macroscopic  polarization/^'  is  modeled  as: 

P‘=R{Elt)  +  Q{El),  (1) 

where  R{El,t)  is  the  time-dependent  reversible  polarization  at  current  time  t  with  P(0,  t)=0 
while  Q{El)  is  the  residual  (irreversible)  polarization.  The  reversible  polarization  is  written  as: 


R*  =  R{El ,  t)  +  £  -^R{El,t-  s)  ^  ds. 


^0  dEl 


where 


R{Elt)  =  R,{El)  +  R,{El) 


1-exp 


ds 

f 

V 


(2) 


(3) 


Both  Rq{EI)  and  Pj(£'3 )  are  function  ofpj .  The  upper  right  superscript  5  indicates  the  history 
time  variable.  The  characteristic  time  r,  measures  the  speed  of  the  polarization  changes  with 
time.  The  irreversible  polarization  is  given  as: 

te‘,  dQ 


s-\: 


dE; 


-dEl 


(4) 


The  rate  of  the  residual  polarization  during  the  polarization  switching  response  is: 
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0, 


Q<El<E^,  dEl  <0  or  -E^<El<  0,  dE^  >  0, 
-E^<El<0,  dEl  <0  or  Q<E[<E^,  dE^  >  0, 


(5) 


dQ‘ 

dEl 


A 


//exp[-ft) 


^3 


], 


-E^<E{<  dEi  <0  or  E<E[<E^,  dEi  >  0, 


V  ^  j 

where  X,  fi,  o),  n  are  the  material  parameters  that  are  calibrated  from  experiments.  E^  and  E^  are 
the  coercive  and  maximum  electric  field,  respectively.  It  is  also  necessary  to  define  a  state  of 
polarization,  which  is  /  ^ ~  the  current  polarization  state. 


analogous  to  the  yield  stress  in  the  overstress  plasticity  theory.  It  is  assumed  that  the  irreversible 
polarization,  Eq.  (9)  is  formed  when  /  )  =  0  and  p'p;  >  0  .  In  a  non-polarized  sample,  the 


current  polarization  state  P^  =  0 ,  and  once  the  ferroelectric  sample  is  completely  polarized  the 
current  polarization  state  is  equal  to  the  saturated  polarization  (P^  =  ox =  -P^ ).  When  the 
materials  undergo  unloading  or  neutral  loading,  E^Pl  <  0 ,  or  loading  that  results  in/(^ 

,  there  is  no  further  update  on  the  polarization  state.  The  other  two  components  of  polarization 
are  expressed  as  P^  =  k^^E[,  and  Pj  =  x-jjp', . 


The  above  constitutive  equations  (Eqs.  (4)-(10))  are  solved  numerically  and  implemented  in 
continuum  based  finite  element.  Discussion  on  the  numerical  algorithm  can  be  found  in  Sohrabi 
and  Muliana  (2013).  The  electro-mechanical  constitutive  model  in  Eqs.  (4)  and  (5)  is  validated 
using  the  experimental  data  of  PZT51  reported  by  Fang  and  Li  (1999).  Figure  10a  shows  the 
polarization  hysteretic  response  (D3-P3)  of  PZT-5 1  subject  to  a  cyclic  electric  field  at  zero  stress. 
The  amplitude  of  the  electric  field  is  1.2  MV/m  with  a  frequency  of  1  Hz.  The  test  started  from 
an  unpolarized  condition  and  with  increasing  the  electric  field,  polarization  takes  place.  The  loop 
in  the  first  cycle  is  higher  by  about  0.03  CIvcX.  After  several  cycles,  the  saturated  polarization 
converges  to  a  constant  value,  which  shows  a  slightly  smaller  value  than  the  one  in  the  first 
cycle.  The  material  parameters  for  the  time-dependent  polarization  in  are  calibrated  from  this 
hysteretic  polarization  response.  The  corresponding  butterfly  hysteretic  responses  during  the  first 
and  saturated  cycles  at  zero  stresses  are  shown  in  Fig.  10b.  It  is  seen  from  Fig.  I  Ob  that  at  the 
coercive  electric  fields  the  strains  from  the  experimental  tests  are  nonzero,  which  is  about  500  ps 
higher.  This  might  be  due  to  the  accumulated  strain  from  the  first  cycle  loading.  In  this  study. 


the  piezoelectric  constant  is  modeled  as  g\j^.  =  gyPP^) 


and  the  parameter  Ci  is 


calibrated  from  the  saturated  butterfly  curve  in  Fig.  11a  after  shifting  the  strain  response  obtained 
from  the  simulation  500  ps  higher.  Tables  2-4  report  the  electro-mechanical  coupling  parameters 
and  elastic  constants,  respectively,  for  PZT-5 1.  The  corresponding  transverse  butterfly  strain 
response  is  shown  in  Fig.  1  lb.  The  nonlinear  electro-mechanical  coupling  model  is  capable  of 
simulating  the  hysteretic  polarization  switching  electro-mechanical  response.  Figures  12  and  13 
shows  the  hysteretic  polarization  and  butterfly  strain  responses  under  cyclic  electric  fields  and 
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constant  compressive  stresses  along  the  electrie  field  lines.  The  model  is  eapable  in  capturing  the 
nonlinear  eleetro-mechanieal  responses  of  PZT. 


Table  2  Material  parameters  for  the  time-dependent  polarization  of  PZT-5 1 


K 

ATO 

K\ 

T\ 

n  p 

CO 

(MV/m) 

(xlO-- 

'F/m) 

(sec) 

(xlO-^F/m) 

(xlO-^F/m) 

0.67* 

70 

225 

1.0 

0.35 

3.0  1.6 

4 

*  from  Fang  and  Li  (1999) 


Table  3  Electro-meehanical  eoupling  parameters  for  PZT-5 1 


^333 

‘^311 

<1 

4 

Pr 

Cl 

m/V) 

(xlO'^F/m) 

(C/m^) 

1520* 

-570* 

38 

42 

0.194* 

0.19 

*  from  Fang  and  Li  (1999) 


Table  4  Elastie  constants  for  PZT-5 1  (Muliana  2010) 


^11  “  ^22 

^33 

Gi2 

(GPa) 

^13  ^23 

Vl2 

Vl3  -V23 

34.48 

33 

13.19 

12.37 

0.307 

0.334 

Figure  10  Hysteresis  polarization  and  butterfly  strain  responses  for  PZT51  at  zero  stress 
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Figure  1 1  Saturated  strain  responses  for  PZT51  at  zero  stress 
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Figure  12  Hysteresis  polarization  responses  under  constant  compressive  stresses  above  the  coercive  stress 
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Figure  13  Butterfly  strain  responses  under  constant  compressive  stresses  above  the  coercive  stress 


The  time-dependent  polarization  model  is  extended  in  order  to  simulate  polarization  switching 
response  of  PZT  5A  at  temperatures:  -75-250°C.  Hooker  (1998)  examined  the  electro¬ 
mechanical  properties  and  polarization  switching  response  of  four  different  PZTs  at  various 
isothermal  temperatures.  Experimental  tests  on  the  polarization  switching  response  were  done  at 
frequency  of  1  Hz  and  at  amplitude  of  electric  field  of  2  MV/m.  In  order  to  determine  material 
parameters  in  the  time-dependent  polarization  model,  for  each  available  data  temperature- 
dependent  material  parameters  were  calibrated  by  matching  the  overall  polarization  switching 
response.  Thus,  it  is  possible  to  have  different  sets  of  material  parameters  that  predict  the 
experimental  data.  Figure  14  illustrates  the  polarization  switching  response  of  PZT-5A  at  various 
isothermal  temperatures.  It  is  seen  that  the  model  is  capable  of  capturing  different  shapes  of  the 
major  hysteresis  curves.  The  calibrated  material  parameters  are  given  in  Table  5,  except  for  the 
values  of  Ec  which  are  obtained  from  Hooker  (1998). 


Table  5  Material  parameters  for  PZT-5A  at  various  temperatures 


T 

(°C) 

Ec 

(MV/m) 

Ko 

(xlO-9 

Ki 

F/m) 

Tl 

(sec) 

X 

(xlO-^  F/m) 

n 

(xlO'^  F/m) 

CO 

-75 

1.41 

10 

10 

1.0 

0.12 

2.0 

0.12 

3.0 

0 

1.2 

20 

25 

1.0 

0.38 

2.0 

0.38 

3.0 

25 

1.15 

20 

35 

1.0 

0.58 

2.0 

0.58 

3.0 

100 

1 

30 

55 

0.5 

1.10 

4.0 

1.10 

5.0 

250 

0.38 

65 

25 

0.1 

2.20 

4.0 

2.20 

6.0 

16 
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Figure  14  Temperature-dependent  polarization  switehing  response  of  PZT-5A  (experimental  data  are 
obtained  from  Hooker  1998) 
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Another  type  of  nonlinear  electro-mechanical  behavior  is  when  the  electric  field  input  is 
prescribed  on  polarized  piezoelectric  materials  with  magnitude  of  electric  field  below  the 
coercive  electric  limit  of  the  materials.  The  model  is  represented  by  the  reversible  polarization 
part  in  Eq.  (7).  A  single  integral  with  a  nonlinear  integrand  is  adopted,  which  is  written  as: 


p  f  £  (•  Se 

£‘  =£(t)=\—(t-s,a(s))a(s)ds+\—(t-s,E(s))E(s)ds 

J  da  •' 


0  ^  0  ^ 

t  t 

D'  =D(0=  f  — f  — (t-5,E(5))E(5)^/5 

da  J  dE 

Alternatively,  one  can  also  consider  the  following  expressions: 

=a(t)=  f— (t-^,£(5))£(5)J5-  f— ^(t-5,E(^))E(^)j5 

J  r)fy  J  5E 


de 

rdD 


=D(t)=  [■^(t-5,£(5))£(5)(i5-l-  f — (t-5,E(s))E(5)(i5 


d£ 


dE 


(11) 

(12) 


(13) 

(14) 


The  time-dependent  kernel  functions  inside  the  above  equations  are  determined  based  on 
experimental  data  on  the  time-dependent  response  of  materials.  In  this  study,  a  series  of 
exponential  function  is  used  for  the  time-dependent  kernel  functions.  A  recursive  integral  model 
is  used  to  solve  the  above  integral  models,  which  is  implemented  at  each  material  point  within 
finite  elements.  Detailed  numerical  procedures  can  be  found  in  Sohrabi  and  Muliana  (2014). 

The  experimental  data  of  a  piezoelectric  ceramics  plate  reported  in  Crawley  and  Anderson 
(1990)  are  first  used  to  validate  the  constitutive  model.  A  cyclic  electric  field  is  applied  through 
the  thickness  direction  of  the  plate  and  the  corresponding  strain  in  the  transverse  direction  is 


measured.  The  strain  from  Eq.  (11)  becomes  =  s) 


ff'jj  dE^(s) 
dE.  ds 


function  is  taken  as  K{t)  =  \- 0.6e  ^ .  Linear  and  nonlinear  functions  are  used  for  the  part. 


ds .  The  time  kernel 


a^ii 


dE, 


ds 

which  are  — ^  =  d 


dE, 


ds 

and  ^  =  d. 


3JJ  ciiiu  -  i*3jj  +  f^^^^^E^E^E^ ,  respectively.  Figure  15  depicts  the  electro- 

3  dE-^ 

mechanical  responses  under  cyclic  electric  field  with  frequency  0.1  Hz  at  three  different 
amplitude  of  electric  fields,  which  are  250,  500,  and  750  V/mm.  Unlike  the  linear  time- 
dependent  electro-mechanical  response,  the  nonlinear  time-dependent  model  can  capture  the 
hysteretic  electro-mechanical  response  of  PZT. 


In  order  to  validate  the  finite  element  implementation  of  the  nonlinear  time-dependent  electro¬ 
mechanical  model  of  polarized  piezoelectric  ceramics,  responses  of  piezoelectric  telescopic 
actuators  are  studied.  Electro-mechanical  responses  of  two  telescopic  actuators  reported  in 
Alexander  et  al.  (2001,  2003)  are  studied.  Figure  16  presents  the  hysteretic  axial  deflection  and 
voltage  input  of  PZT-568  actuator.  Frequency  of  O.lHz  is  considered  in  the  simulation.  Slight 
nonlinear  response  is  shown.  The  electric  potential  and  displacement  fields  at  maximum  and 
minimum  voltages  in  the  actuator  is  also  shown.  Another  actuator,  MSI-53,  is  considered,  whose 
electro-mechanical  response  under  electric  field  input  is  illustrated  in  Fig.  17.  The  response 
shows  mainly  linear  hysteretic  behavior,  which  is  most  likely  due  to  relatively  small  voltage 
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applied.  The  contours  of  axial  displacement  and  electric  potential  are  also  given  in  Fig.  17. 
Overall  the  time-dependent  electro-mechanical  model  is  capable  in  capturing  the  hysteretic 
response  of  piezoelectric  materials  and  structures. 


Figure  15  Linear  (left)  and  nonlinear  (right)  time-dependent  electro-mechanical  response  of  PZT  under 
cyclic  electric  field. 


Figure  16  Electro-mechanical  response  of  PZT-568  telescopic  actuator 


■  300[V] 

Electric  Potential  ® 

■  -300  [V] 


20  [pm] 


Axial  Deformation 


10  [p  m] 
0[pm] 


Emax=300  [V] 


-20  [pm] 

-10  [p  m] 
0[pm] 


Emin=-300  [V] 


Figure  17  Electro-mechanical  response  of  MSI-53  telescopic  actuator 
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In  order  to  be  able  to  reasonably  predict  the  overall  behaviors  of  smart  structures  under  various 
thermo-electro-mechanical  stimuli,  it  is  necessary  to  incorporate  responses  of  the  constituents  in 
the  smart  composites.  Two  main  constituents,  which  are  polymers  and  piezoelectric  ceramics, 
are  considered  for  the  smart  composites.  Since  polymers  and  piezoelectric  ceramics  often 
experience  nonlinear  and  time-dependent  responses  under  various  histories  of  thermo-electro- 
mechanical  external  stimuli,  nonlinear  and  time-dependent  constitutive  models  are  developed  for 
polymers  and  piezoelectric  ceramics.  The  nonlinearity  is  considered  due  to  high  magnitude  of 
stimuli.  Numerical  methods  are  also  developed  for  obtaining  solutions  to  the  nonlinear  and  time- 
dependent  constitutive  models,  which  are  implemented  within  finite  elements.  The  constitutive 
models  are  also  integrated  to  several  micromechanical  models  in  order  to  study  performance  of 
smart  composites.  For  convenience  in  solving  the  nonlinear  time-dependent  constitutive 
equations,  an  incremental  formulation  is  used  and  within  each  time  increment  a  linearized 
relation  is  considered.  The  linearized  relation  for  a  thermo-electro-mechanical  response  is: 

z  =  nr 


r'”'  =  |A£-jj,A£'22,  A£-33,A£-j2,A£'j2,A£'23,AE'j,AE'2,AE'3,Ar} 

Z'*'  =  {Acrjj,Acr22,Acr33,Acrj2,  Af7j2,Acr23,ADj,  ADj,  ADj,  AA} 


C 


a  = 


e 


P 


-h 

g 

P 


(15) 


where  Q  is  the  consistent  tangent  matrix  that  relate  the  incremental  field  variables.  In  case  of 
nonlinear  response,  an  iterative  method  is  also  included  within  the  incremental  formulation. 


4.3)  Micromechanical  Models 

Several  micromechanical  models  for  fiber  and  particulate  composites  are  formulated  in  order  to 
determine  the  overall  (macroscopic)  response  of  composites  due  to  conduction  of  heat,  thermo¬ 
mechanical  inputs,  and  thermo-electro-mechanical  stimuli.  The  constitutive  models  discussed  in 
Section  4.2  are  implemented  for  the  constituents  in  the  micromechanical  models.  First, 
simplified  unit-cell  models  are  formulated  for  fiber  and  particle  composites,  shown  in  Fig.  18a 
and  b.  The  simplified  unit-cell  models  consist  of  several  subcells.  In  case  of  fiber  composites, 
four  fiber  and  matrix  subcells  are  considered,  while  for  particle  composites,  there  are  eight 
subcells.  Other  idealized  micromechanical  models  include  detailed  cylindrical  fiber  of  circular 
cross-section  and  solid  spherical  particle  arrangements  in  homogeneous  matrix.  Uniform  and 
random  distributions  of  the  inclusions  (fibers  and  particles)  are  considered.  The  overall  responses 
determined  from  the  unit-cell  models  are  compared  to  the  ones  obtained  from  the 
micromechanical  models  with  detailed  fiber  and  particle  arrangements,  and  also  with  available 
experimental  data.  The  unit-cell  model  for  particle  reinforced  composites  is  also  used  to  study 
overall  performance  of  functionally  graded  materials.  The  third  micromechanical  model  for 
particle  composites  incorporate  detailed  microstructural  morphologies,  i.e.,  shapes  of  inclusions, 
pores,  and  distributions  of  inclusions  in  the  matrix,  shown  in  Fig.  18c.  This  is  done  to  mimic 
more  realistic  microstructures  and  examine  the  effect  of  detailed  microstructural  morphologies  of 
composites  on  the  overall  macroscopic  responses  of  composites.  In  this  project,  perfectly  bonded 
interfaces  between  inclusions  and  matrix  are  considered. 
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a)  Fiber-reinforced  composite 


Idealized  microstructure  Idealized  microstructure 

c)  Micromechanics  models  based  on  SEM  microstructural  images 


Figure  18  Several  micromechanical  models  for  composites 


4.3.1)  Heat  conduction  and  deformation  analyses 

Micromechanical  models  for  fiber  and  particle  composites  are  formulated  in  order  to  determine 
the  overall  response  of  composites  undergoing  heat  conduction  and  deformations.  Sequentially 
and  fiilly  coupled  heat  conduction  and  deformation  analyses  are  considered  at  multiple  length 
scales.  In  the  sequentially  coupled  analysis  it  is  assumed  that  the  temperature  field  influences  the 
deformation  field;  thus,  the  governing  equations  for  the  conduction  of  heat  are  solved  first  in 
order  to  determine  the  temperature  field  in  the  composite  bodies  followed  by  solving  the 
governing  equations  for  the  deformation.  In  the  case  of  fiilly  coupled  thermo-mechanical 
analysis,  the  temperature  field  affects  the  deformations  in  the  composite  bodies  and  the 
dissipation  of  energy  from  the  viscoelastic  responses  also  influence  the  temperature  field.  In  such 
situation,  coupled  governing  equations  for  the  heat  conduction  and  deformation  are  solved 
concurrently  at  multiple  length  scales.  It  is  noted  that  the  problems  considered  here  are  under 
small  deformation  gradients,  in  which  a  linearized  strain  measure  is  considered,  and  the  thermal 
and  mechanical  properties  of  the  constituents  are  allowed  to  vary  with  temperatures,  stresses,  and 
time,  leading  to  nonlinear  responses  of  the  composites. 
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The  unit-cell  models  for  fiber  and  particle  reinforced  composites  are  first  considered  in  order  to 
determine  the  macroscopic  thermo-mechanical  behaviors  of  homogenized  composites.  The  unit¬ 
cell  models  consist  of  inclusion  and  matrix  constituents,  which  are  represented  by  subcells,  and 
the  nonlinear  thermo-mechanical  constitutive  models  discussed  in  Section  4.2  are  used  for  the 
subcells.  The  unit-cell  models  are  implemented  in  continuum  finite  elements  and  used  to  study 
temperature  and  deformation  fields  in  the  homogenized  composite  medium.  The  second  type 
micromechanical  model  treats  heterogeneous  composites  explicitly.  Heterogeneous  composites 
that  consist  of  inclusions  (particles  or  fibers)  randomly  distributed  in  homogeneous  matrix  are 
generated  using  3D  continuum  elements  in  FE  framework.  For  each  volume  fraction  (VF)  of 
inclusions,  the  models  of  heterogeneous  composites  with  different  sizes  of  inclusions  and 
arrangements  are  generated  such  that  these  models  represent  more  realistic  volume  elements. 


In  deriving  the  unit-cell  models,  the  following  assumptions  are  made.  The  microstructural  scale, 
which  is  represented  by  the  unit-cells,  is  much  smaller  than  the  macroscopic  scale  so  that  the 
transient  heat  conduction  at  the  microstructural  scale  can  be  neglected.  Each  inclusion  is  fiilly 
surrounded  by  homogeneous  matrix  and  perfectly  bonded  to  the  matrix.  In  each  subcell  within 
the  unit-cell,  the  variations  of  field  variables  are  assumed  uniform.  The  effective  field  variables 
for  the  unit-cell  models  are  determined  based  on  the  volume  averaging  scheme  (Aboudi  2005): 

q  =  -^  I  (16-1) 

^  m=\  vim)  ^ 


^  m=\  yim)  V  m=\ 


(16-2) 


(16-3) 


m=l  T/(m) 


lyy 


(m)g(m) 


(16-4) 


where  q,  cp,  a,  and  8,  are  the  heat  flux,  temperature  gradient,  stress,  and  strain  fields, 
respectively.  An  overbar  denotes  the  average  (macroscopic)  response  of  the  unit-cell, 
denotes  the  volume  fraction  of  each  subcell,  the  total  volume  fraction  of  the  unit-cell  is 


E  =  ^  and  N  is  the  number  of  sucells.  Equations  (16-l)-(16-4)  are  exact  only  for  linear 


thermo-elastic  responses  of  composites.  Due  to  the  nonlinear  and  time-dependent  responses  in 
the  constituents,  the  solutions  for  the  deformation  fields  are  determined  incrementally.  The 
incremental  forms  of  the  effective  field  variables  at  the  current  time  are  q'  =  q'"^'  -i-  Aq' , 

(p'  =  (p'"^'  +  Aip' ,  ct'  =  +  Aa‘  and  T'  =  +  As' .  The  superscript  (t-At)  denotes  previous 

converged  solutions  and  the  prefix  A  indicates  the  incremental  variables  within  a  current  time 
increment  At.  A  linearized  relation  is  first  imposed  to  the  incremental  parts  of  the  field  variables. 
An  iterative  corrector  scheme  is  then  performed  in  order  to  satisfy  both  the  micromechanical 
constraints  and  the  nonlinear  constitutive  equations. 


The  micromechanical  models  and  constitutive  equations  are  implemented  within  finite  elements, 
in  which  displacement  and  temperature  fields  are  defined  at  the  nodes.  Using  the  spatial 
derivatives  of  the  shape  functions,  the  average  temperature  gradient  and  strain  are  determined  at 
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each  material  (Gaussian)  point  within  finite  elements.  The  macroscopic  (average)  temperature 
gradient  and  strain  are  related  to  the  temperature  gradient  and  strain  in  each  subcell  in  the  unit¬ 
cell  models  by  the  following  relations: 

AtpW  =  and  (17) 

where  the  second  order  concentration  tensor  relates  the  temperature  gradient  in  each  subcell 

to  the  average  temperature  gradient  of  the  unit-cell  and  the  fourth  order  concentration  tensor 

relates  the  strain  in  each  subcell  to  the  average  strain  of  the  unit-cell.  In  order  to  define  the 
concentration  tensors,  micromechanical  relations  and  constitutive  equations  for  each  constituents 
are  used,  which  are  discussed  in  Muliana  and  Kim  (2010),  Muliana  and  Lin  (201 1),  Khan  and 
Muliana  (2012b).  Using  linearized  constitutive  relations  for  the  heat  conduction  and  mechanical 
deformation,  the  relations  between  the  micro-macro  scale  field  variables,  and  the  volume 
averaging  scheme,  the  effective  linearized  constitutive  relations  are: 

Aq  «  —  Z  (18-1) 

^  m=\ 

Ao  «  —  Z  (18-2) 

^  m=\ 

and  the  effective  linearized  thermal  conductivity  and  stiffness  are  k  «  —  Z  and 

V  tn=\ 

C  «  —  Z  ,  respectively.  It  is  noted  that  when  nonlinear  behaviors  are  considered, 

^  m=\ 

the  linearized  micromechanical  models  are  used  as  trial  (predictor)  solutions,  and  an  iteration  is 
added  in  order  to  minimize  residual  from  linearizing  the  nonlinear  responses. 

The  following  cases  have  been  considered  in  order  to  examine  the  average  (macroscopic) 
response  of  the  homogenized  composites  obtained  from  the  unit-cell  models. 

a)  Heat  conduction  in  fiber  reinforced  composites 

The  unit-cell  model  with  four  subcells  that  represent  fiber  reinforced  polymer  (FRP)  composites 
is  used  to  analyze  heat  conduction  in  composites.  Experimental  data  available  in  literature  are 
used  to  validate  the  effective  thermal  conductivity  of  fiber  composites  (see  Muliana  and  Kim 
2010  for  a  detailed  discussion).  Figure  19  depicts  the  effective  thermal  conductivities  in  the  axial 
and  transverse  the  fiber  directions  of  carbon/epoxy  composites.  Experimental  data  are  obtained 
from  Mclvor  et  al.  (1990).  Good  predictions  are  observed. 

Next,  heat  conduction  behaviors  obtained  from  the  unit-cell  model  are  compared  to  the  ones  of 
heterogeneous  bodies  incorporating  detailed  fiber  arrangements.  Figure  20  illustrates  the  finite 
element  meshes  for  homogenized  FRP  and  heterogeneous  FRP  with  detailed  fiber  arrangements. 
In  case  of  heterogeneous  FRP,  random  and  uniform  fiber  arrangements  are  considered.  Figure  21 
shows  the  temperature  fields  during  heat  conduction  in  the  axial  fiber  direction  for  composites 
with  50%  fiber  volume  content.  Relatively  small  fiber  cross-sectional  area,  i.e.,  0.1  mm^,  is 
considered  in  the  heterogeneous  FRP.  Figure  22  shows  the  comparisons  between  temperature 
fields  in  the  homogenized  and  heterogeneous  FRP  composites  when  a  larger  fiber  cross-sectional 
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area  is  considered.  The  steady  state  time  for  the  transient  heat  conduction  in  the  axial  and 
transverse  fiber  directions  determined  from  the  homogenized  and  heterogeneous  composites  with 
several  fiber  volume  fractions  are  summarized  in  Table  6.  The  heterogeneous  response  is 
denoted  by  microstructural  fiber  (MF)  model.  The  central  processing  unit  (CPU)  times  needed  to 
complete  the  heat  conduction  analyses  to  reach  steady  state  for  both  models  are  also  reported.  As 
expected,  the  unit-cell  model  requires  significantly  less  CPU  time  to  complete  the  analyses. 


Teniperatme  (“q  Tsmperature  (”q 

Figure  19.  Temperature-dependent  thermal  conductivity  of  GY80  carbon/69epoxy 


FE  mesh  of  heterogeneous  FE  mesh  of  homogenized  FRP 


model 

Figure  20.  FE  microstructures  for  heterogeneous  and  homogenized  FRP  systems 
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Figure  21  Temperature  profile  due  to  eonduetion  along  the  fiber  direetion 
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Figure  22  Temperature  profile  due  to  conduction  along  the  fiber  direction 
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Table  6  Times  to  steady  state  for  transient  heat  transfer  analyses,  eomparison  of  the  two  models 


Direction  of  heat 
conduction 

Model 

Fiber  Volume  Fraction 

10 

20 

40 

50 

67 

75 

Steady  State  time 
(Longitudinal)  [s] 

UC  model 

220 

191.5 

157.5 

146 

131.5 

126 

CPU  time  [s] 

781.5 

688.1 

565.8 

526.5 

473.8 

459.9 

MF  model 

219 

189 

153.5 

143.5 

129 

124 

CPU  time  [s] 

8686 

15833 

8666 

4150 

5175 

4987 

%  Difference  to  reach  steady  state 

0.46% 

1.32% 

2.61% 

1.74% 

1.94% 

1.61% 

Steady  State  time 
(Trans)  [s] 

UC  model 

263 

253.5 

226.5 

209.5 

179 

163 

CPU  time  [s] 

935.5 

911.6 

810.5 

754.1 

642.5 

590.3 

MF  model 

254 

240 

214 

206.5 

173.5 

154 

CPU  time  [s] 

10077 

20175 

12045 

5968 

6916 

6091 

%  Difference  to  reach  steady  state 

3.54% 

5.63% 

5.84% 

1.45% 

3.17% 

5.84% 

b)  Coupled  heat  conduction  and  deformation  in  particle  reinforced  composites 

In  this  section  sequentially  and  fully  coupled  heat  conduction  and  deformation  are  presented  for 
particle  reinforced  composites.  The  unit-cell  model  is  used  to  determine  the  average 
(macroscopic)  response  in  the  case  of  homogenized  composites  are  considered.  The  responses 
are  also  compared  to  the  heterogeneous  model  with  detailed  particle  arrangements.  Figure  23 
shows  the  homogenized  and  heterogeneous  composites.  In  the  case  of  heterogeneous  model  is 
considered,  different  particle  sizes  are  also  considered  for  composites  at  fixed  particle  volume 
contents.  The  composites  are  subjected  to  temperature  gradient  on  the  opposite  left  and  right 
surfaces  and  the  left  surface  with  lower  temperature  is  constraint  so  that  the  displacement  is  zero. 
Figures  24  and  25  show  the  temperature  and  displacement  profiles  during  the  transient  heat 
conduction  process.  Responses  from  the  unit-cell  model  (homogenized  composites)  are 
compared  to  the  ones  of  the  heterogeneous  bodies.  The  corresponding  thermal  stresses  during 
heat  conduction  are  also  presented  in  Fig.  26.  As  expected  the  homogenized  composites  give 
smooth  variations  in  thermal  stresses,  while  the  heterogeneous  bodies  can  capture  stress 
discontinuities  and  variations  at  the  interfaces  between  particles  and  matrix.  The  unit-cell  model 
is  also  used  to  predict  fully  coupled  heat  conduction  and  deformation  in  composites  and 
functionally  graded  systems,  as  depicted  in  Fig.  7. 


(a) 


(b) 


Figure  23  Finite  element  meshes  of  a)  homogenized  and  b)  heterogeneous  particle  composites. 
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Figure  24  Temperature  profiles  for  eomposites  with  30%  partiele  eontent:  a,  b,  and  e  are  aetual  values 
measured  at  the  top  surfaee;  d,  e,  and  f  are  mean  values  measured  from  top  and  bottom  surfaees 
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Figure  25  Displacement  profiles  for  composites  with  30%  particle  content:  a,  b,  and  c  are  actual  values 
measured  at  the  top  surface;  d,  e,  and  f  are  mean  values  measured  from  top  and  bottom  surfaces 
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Figure  26  Thermal  stress  in  the  direction  of  temperature  gradient  during  the  transient  process  for 
composites  with  30%  particle  content  and  different  particle  sizes. 


c)  Viscoelastic  and  viscoplastic  responses  of  particle  and  fiber  reinforced  composites 

The  unit-cell  models  for  the  particle  and  fiber  composites  are  now  used  to  study  the  viscoelastic 
and  viscoplastic  responses  of  composites.  The  constitutive  models  for  viscoelastic-viscoplastic 
responses  are  integrated  to  unit-cell  models  and  also  to  finite  elements.  The  overall  responses 
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from  the  unit-cell  models  are  compared  to  the  ones  determined  from  the  composites  with  detailed 
arrangements  of  inclusions  and/or  with  experimental  data  available  in  literature. 

Figure  27  shows  the  creep-recovery  behaviors  of  particle  composites  from  the  unit-cell 
(homogenized  composite)  and  composite  with  detailed  particles.  The  viscoplastic  models  of 
Valanis  and  Perzyna  are  also  considered.  It  is  seen  that  the  unit-cell  model  can  capture  the 
average  (macroscopic)  viscoelastic-viscoplastic  responses  of  the  composites.  The 
micromechanics  models  are  also  use  to  study  a  uniaxial  tensile  response  of  a  composite  bar.  Both 
unit-cell  model  and  microstructural  model  with  detailed  particle  arrangements  are  considered.  As 
seen  in  Fig.  28,  the  homogenized  responses  can  predict  the  overall  tensile  behavior  of  the 
composites.  The  stress  (von  Mises)  fields  of  the  two  composite  models  under  a  uniaxial  tensile 
loading  at  different  times  are  illustrated  in  Fig.  29.  The  limitation  in  the  homogenized  model  is  in 
capturing  detailed  stress  variations  and  concentrations  in  the  composites  as  the  unit-cell  model 
gives  the  average  quantities  of  field  variables. 
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Figure  27  Creep  reeovery  responses  of  partiele  eomposites  with  different  partiele  volume  fraetion 


Figure  28  Creep  responses  under  uniaxial  tension  of  30  MPa 
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a)  Homogenized  model 


t=0.0001  sec 


b)  Heterogeneous  model 


t=0.0001  sec 


t=2070 

H  60  MPa  □  30  MPa  H  0 

Figure  29  von-Misses  stress  distribution  of  creep-recovery  responses  under  30  MPa 
a)  homogenized  model  b)  heterogeneous  model 

Next,  the  unit-cell  model  for  fiber  composites  is  used  to  study  viscoelastic-viscoplastic  responses 
of  composites.  Experimental  data  on  the  viscoelastic-viscoplastic  response  of  glass  FRP 
unidirectional  composites  reported  by  Megnis  and  Varna  (2003)  are  used  to  validate  the  unit-cell 
model.  The  tests  were  performed  on  several  off-axis  composite  specimens:  20°,  30°,  45°,  90°  and 
0°  and  stress  levels.  The  viscoplastic  parameters  in  the  Perzyna’s  model  were  characterized  by 
fitting  the  shear  strain  experimental  data  of  the  off-axis  specimens  at  high  stress  levels.  The 
remaining  data  of  the  off-axis  specimens  that  were  not  used  in  the  calibration  are  used  to  validate 
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the  unit-cell  model.  Figure  30  shows  the  strain  responses  along  the  loading  and  transverse  to  the 
loading  directions  for  composites  with  90°  off-axis  fiber  angle.  Responses  from  the  90° 
composite  show  mainly  viscoelastic  behaviors.  Figure  31  presents  the  responses  for  20°  off-axis 
fiber  and  significant  viscoplastic  response  are  seen  in  the  shear  component.  Detailed  discussion 
on  the  model  and  material  calibration  can  be  found  in  Jeon  and  Muliana  (2012). 

The  responses  from  the  homogenized  composite  model  are  also  compared  to  those  from  the 
micromechanics  models  with  detailed  fiber  arrangements  (heterogeneous  model),  as  illustrated  in 
Fig.  20.  Uniform  and  random  fiber  arrangements  are  considered  in  the  heterogeneous  model.  For 
this  simulation,  carbon  fiber  is  used  for  the  fiber  constituent  and  the  high-density  polyethylene 
(HDPE)  is  used  for  the  matrix.  Creep  recovery  behaviors  of  homogenized  and  heterogeneous 
composites  with  different  fiber  volume  fractions  are  shown  in  Fig.  32.  Loading  is  prescribed  in 
the  transverse  fiber  direction,  and  the  corresponding  strain  outputs  in  the  transverse  and  axial 
directions  are  compared.  Again,  the  unit-cell  model  which  represents  the  homogenized  (average) 
behaviors  of  the  composites  can  capture  the  overall  time-dependent  and  inelastic  responses. 
Tables  7  and  8  summarize  percent  differences  in  the  strain  responses  calculated  from  the  unit¬ 
cell  and  micromechanical  model  with  detailed  fiber  arrangements,  which  in  this  study  the  largest 
difference  is  around  7%.  One  of  the  main  advantages  of  using  the  unit-cell  model  is  in  reducing 
the  computational  cost  since  the  stress  and  strain  fields  are  sampling  at  several  subcells  instead 
of  over  the  entire  domain  of  composites.  Table  9  presents  the  central  processing  unit  (CPU)  time 
needed  in  analyzing  the  time-dependent  responses  of  composites  using  the  unit-cell  and  detailed 
micromechanical  models. 
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Figure  30  Creep  strain  for  [90] m  laminate.  (  o)  experiment  strain  data,  ( — )ealeulated  viseoelastie  strain. 

Applied  stress=50MPa  (Transverse  of  fiber  direetion). 
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Figure  31  Creep  strain  for  [20]  m  laminate. 


Table  7  Percent  difference  of  the  strains  from  the  unit-cell  and  FE  micromechanical  models 


Fiber  Array 

Time  [s] 

Loading 

Direction 

Strain 

Difference  % 

[10%  Vf] 

[20%  Vf] 

[50%  Vf] 

Unit-cell  Model  / 
Detailed  FE  Model 

VE 

1800 

Fiber 

Sll 

2.43% 

2.49% 

4.45% 

£22 

2.05% 

1.82% 

5.52% 

Transverse 

£ll 

1.67% 

1.79% 

4.86% 

£22 

1.99% 

2.03% 

5.37% 

VE+VP 

1800 

Fiber 

£ll 

2.42% 

2.49% 

4.46% 

£22 

2.05% 

1.82% 

5.54% 

Transverse 

£ll 

0.02% 

6.31% 

5.02% 

£22 

0.31% 

3.90% 

5.69% 
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Table  8  Percent  difference  of  the  strains  from  the  unit-cell  and  FE  micromechanical  models  with 
uniform  detailed  fiber  arrangements 


Loading  Direction 

Model 

Vf40 

Vf54 

Vf  67 

Vf  75 

Axial  Loading 

sii  [mm/mm] 
(1800  s) 

Unit-Cell 

L08E-04 

8.01E-05 

6.51E-05 

5.79E-05 

Detailed  FE 

LllE-04 

8.21E-05 

6.68E-05 

5.94E-05 

%  Difference 

2.53% 

2.54% 

2.53% 

2.53% 

Transverse 

Loading 

822  [mm/mm] 
(1800  s) 

Unit-Cell 

L73E-03 

1.26E-03 

9.17E-04 

7.27E-04 

Detailed  FE 

L77E-03 

1.34E-03 

9.87E-04 

7.58E-04 

%  Difference 

2.47% 

6.05% 

7.12% 

4.10% 

Number  of  Elements 

244590 

240630 

234360 

247290 

Fiber  Arrangement 

12*12 

13*15 

15*16 

15*18 

Arrangement  Ratio 

1.000 

0.867 

0.938 

0.833 

Table  9  CPU  time  for  VE  and  VE+VP  in  homogeneous  and  heterogeneous  models 


FE  Models 

Vf  [%] 

Total 

Elements 

CPU  Time  [s] 

Fiber  Dir. 

Transverse  Dir. 

VE 

Unit-cell  model 

10 

1 

2.1 

75.4 

20 

1 

2.0 

75.6 

50 

1 

1.6 

168.7 

Detailed  FE  model 
[Uniform  Array] 

10 

73728 

41170  [11  Hr] 

534776  [149  Hr] 

20 

73728 

10952  [3  Hr] 

164220  [46  Hr] 

50 

137490 

25724  [7  Hr] 

20632  [6Hr] 

VE+VP 

Unit-cell  model 

10 

1 

1.9 

78.6 

20 

1 

1.9 

79.7 

50 

1 

1.8 

3836.6 

Detailed  FE  model 
[Uniform  Array] 

10 

73728 

92753  [26  Hr] 

11059224  [128  Days] 

20 

73728 

137217  [38  Hr] 

2095620  [24  Days] 

50 

137490 

24582  [7  Hr] 

3539790[41  Days] 
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Figure  32  Creep  recovery  under  loading  in  the  transverse  fiber  direction.  Top:  response  is  the  loading 
direction  and  Bottom:  response  in  the  perpendicular  direction  to  loading 


The  unit-cell  model  is  also  used  to  study  the  thermo-mechanical  responses  of  FRP  composites. 
Examples  of  prediction  of  the  coefficients  of  thermal  expansion  in  the  axial  and  transverse  fiber 
directions  are  shown  in  Fig.  33.  Experimental  data  are  obtained  from  Karadeniz  and  Kumlutas 
(2007).  The  unit-cell  model  can  capture  the  thermal  expansion  in  the  axial  fiber  direction,  while 
some  mismatches  are  seen  in  the  transverse  fiber  directions.  The  reason  of  the  mismatches  are 
due  to  the  thermal  stress  effects  in  the  constituents  due  to  temperature  changes.  The  unit-cell 
model  is  limited  in  capturing  the  detailed  stress  distributions  in  the  constituents. 

Figure  34  illustrates  the  unit-cell  predictions  of  stress  relaxation  behaviors  in  AFCs  at  different 
temperatures  and  strain  levels.  The  AFCs  consist  of  piezoelectric  ceramic  long  fibers  and  epoxy 
matrix.  The  relaxation  responses  are  due  to  the  viscoelastic  matrix,  which  are  more  pronounced 
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at  elevated  temperatures  and  higher  strain  levels.  The  unit-cell  model  is  capable  in  capturing  the 
relaxation  responses  of  AFCs  when  damage/degradation  is  absent.  A  significant  mismatch 
between  the  model  and  experimental  data  is  observed  under  a  strain  of  1.61%,  in  which 
damage/cracking  in  the  brittle  ceramic  fiber  is  observed  in  the  specimen. 


Figure  33  Thermal  expansion  coefficients  in  the  axial  (left)  and  transverse  (right)  fiber  directions. 


Figure  34  Stress  relaxation  responses  of  AFC  at  different  temperatures  and  strain  levels. 


4.3.2)  Thermo-electro-mechanical  analyses 

Micromechanical  models  for  fiber  and  particle  composites  are  formulated  in  order  to  determine 
the  overall  thermo-electro-mechanical  response  of  composites  undergoing  heat  conduction  and 
electro-mechanical  stimuli.  Unit-cell  models  are  first  considered.  A  general  linearized 
constitutive  relations  written  in  Eq.  (15)  is  adopted  and  used  within  a  volume  averaging  scheme 
in  the  unit-cell  model.  Micromechanical  models  with  detailed  fiber  and  particle  arrangements  are 
also  generated  and  used  to  study  the  thermo-electro-mechanical  response  of  active  composites. 
Furthermore,  the  responses  from  the  unit-cell  models  are  compared  to  currently  available 
micromechanics  models  such  as  Mori-Tanaka  and  Self  Consistent  models. 


Following  a  linearized  constitutive  model  in  Eq.  (15),  the  average  (effective)  field  variables, 
indicated  by  an  over-bar,  in  the  unit-cell  models  are  defined  as: 
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(19-1) 


r  =  -^  j  «— 

^  m=\  y{m)  V 

Z  =  — ^  I  Z^"‘\x)dV^'"^ 

^  ^  ^  m=\ 


(19-2) 

The  variables  and  are  the  eomponents  of  the  average  field  quantities  in  eaeh  subeell. 


J 

^  m=\  y(m) 


N 

The  unit-cell  volume  V  is  defined  by  V  =  .  This  study  deals  with  nonlinear  electro- 

m=l 

mechanical  response  for  the  inclusions  and  viscoelastic  response  of  the  matrix  constituent,  which 
requires  incorporating  history  of  loading  to  determine  the  overall  response  of  the  composites. 

The  solutions  for  the  field  variables  are  performed  incrementally.  The  incremental  forms  of  the 
effective  and  local  field  variables  at  current  time  are: 


r' 


— t-At  — t 

=  r  +Ar 


— t  — t-At  — t 


z 

r 

z 


=  Z  +AZ 

(m),t  _  j-n(m),f-Af 


(20) 


where  the  superscript  t-At  denotes  the  previous  time  and  all  the  field  variables  associated  with  t- 
At  are  considered  as  history  variables.  At  each  time  increment  At,  it  is  necessary  to  calculate  the 
incremental  field  variables  at  current  time.  A  linearized  micromechanical  relation  is  imposed  to 
the  incremental  parts  of  the  field  variables  and  a  residual  arouse  from  the  linearization  is 
corrected  incrementally  to  satisfy  both  micromechanical  constraints  and  constitutive  equations 
for  the  inclusions  and  matrix.  For  simplicity,  the  superscript  related  to  time  is  dropped  and  the 
rest  of  the  incremental  formulation  are  presented  at  current  time  t.  The  local  average  independent 
field  variables  can  be  expressed  in  terms  of  effective  (homogenized)  field  variables  which  in  the 
incremental  form  are  written  as: 

ArC”)  =  F(”’)Af  (21) 

where  is  the  concentration  tensor  for  the  subcell  (m). 


The  above  averaging  schemes  give  an  exact  effective  response  of  the  homogenized  composites 
only  when  linear  constitutive  equations  are  used  for  all  constituents  (subcells)  in  the  unit-cell 
model.  In  the  case  of  nonlinear  constitutive  models  are  considered  for  the  constituents,  imposing 
the  linearized  micromechanical  relations  for  the  incremental  field  variables  could  lead  to 
residual.  At  every  instant  of  time  a  corrector  scheme  is  added  to  minimize  the  residual.  This 
study  uses  the  Newton-Raphson  iterative  method  as  a  corrector  scheme.  Once  the  residual  has 
been  minimized,  the  components  of  the  concentration  matrices  are  determined  and  the  field 

variables  Ar^™*  in  each  subcell  due  to  the  prescribed  boundary  conditions  at  the  macro 
(homogenized  composite)  level  are  calculated  using  Eq.  (21).  Finally,  the  constitutive  equations 
in  Section  4.2  are  used  to  obtain  the  field  variables  AZ^™*  in  each  subcell  and  the  effective  field 
variables  at  every  instant  of  time  are  evaluated  using  the  volume  average  scheme  in  Eq.  (19). 
From  the  micro-macro  scale  relation  in  Eq.  (21)  and  constitutive  relation  for  each  subcell  in  Eq. 
(15),  the  effective  field  variables  from  Eq.  (19)  is  given  as: 
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(22) 


Az«  —  y  p^('«)n('«)F(™)Ar 

_  Y  N 

and  the  effective  linearized  properties  are  H  ^ 

^  m=\ 


a)  Nonlinear  response  of  polarized  active  inclusions 

First  the  effective  responses  of  active  fiber  composites  under  a  uniaxial  tensile  creep  of  10  MPa 
are  simulated  using  the  unit-cell  model.  Micromechanical  models  with  detailed  fiber 
arrangements,  which  are  implemented  in  FE,  are  also  considered  for  comparisons.  Figures  35 
and  36  illustrate  the  corresponding  axial  and  transverse  strains  and  electric  flux  in  PZT- 
7A/LaRC-Si  composites  during  creep  loading  at  different  temperatures.  The  viscoelastic  matrix 
induces  time-dependent  changes  in  the  electric  flux  of  the  composites.  The  unit-cell  model  is 
capable  in  capturing  the  overall  thermo-electro-mechanical  responses  of  active  composites. 

Next,  field  variables,  i.e.,  temperature,  displacement,  stress,  electric  field,  and  electric  flux,  in 
active  fiber  during  transient  heat  conduction  throughout  the  composite  body  are  examined  using 
the  unit-cell  model.  The  unit-cell  model  allows  incorporating  different  time-  and  field  dependent 
behaviors  for  the  constituents  during  the  conduction  of  heat  and  deformation  in  the  composite 
body.  Consider  an  active  fiber  composite  panel  where  the  longitudinal  fibers  are  aligned  in  the  X3 
direction  undergoing  conduction  of  heat  along  the  longitudinal  fiber  axis.  It  is  assumed  that  the 
composite  panel  is  perfectly  insulated  so  that  the  amount  of  heat  loss  is  ignored.  The  analytical 
solutions  are  obtained  by  imposing  a  plane  stress  condition.  Figure  37  depicts  the  axial 
displacements,  stresses,  electric  fluxes,  and  electric  fields  at  several  instants  of  time  during  the 
transient  heat  conduction.  The  non-uniform  temperature  profiles  at  early  times  lead  to  relatively 
high  axial  stress,  electric  field,  and  electric  flux.  As  time  progresses  and  steady  state  condition  is 
reached,  the  axial  stress  and  electric  field  should  drop  to  zero,  while  the  electric  flux  approaches 
a  constant  value.  On  the  other  hand,  the  transverse  stress  increases  as  time  progresses  and  the 
maximum  value  is  at  the  steady  state  condition,  which  is  expected. 

Next,  the  unit-cell  models  for  electro-active  composites  are  compared  to  available  experimental 
data  and  other  micromechanics  models,  i.e.,  Mori  Tanaka  (MT)  and  Self-consistent  (SC)  models. 
Detailed  formulations  of  the  MT  and  SC  models  with  nonlinear  electro-mechanical  constitutive 
models  for  the  inclusions  are  given  in  Lin  and  Muliana  (2013).  The  overall  electro-mechanical 
properties  determined  using  the  micromechanics  models  are  first  compared  with  existing 
experimental  data  with  linear  piezoelectric  properties.  Chan  and  Unsworth  (1989)  experimentally 
examined  a  piezoelectric  FRC  with  continuous  fiber  reinforcement  PZT-7A  embedded  in  an 
Araldite  D  medium.  Figure  38  illustrates  the  comparisons  of  effective  electrical  and  mechanical 
properties  of  active  fiber  composites.  All  micromechanics  models  are  in  agreements  with  the 
experimental  data  over  the  whole  range  of  volume  fractions.  Figure  39  depicts  the 
micromechanics  models  predictions  of  the  piezoelectric  particle  composites,  which  was 
experimentally  studied  by  Furukawa  et  al.  (1976).  The  piezoelectric  composite  consists  of 
spherical  PZT-5  particles  embedded  in  an  epoxy  medium.  Finally,  Fig.  40  shows  the  linear  and 
nonlinear  electro-mechanical  responses  determined  from  the  UC  and  MT  models.  As  seen  in  the 
figure,  at  high  electric  field  input,  the  nonlinear  responses  lead  to  higher  strains  in  the 
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Axial  strain  (833%)  Axial  strain  (833%)  Axial  strain  (633%) 


composites.  The  micromechanics  models  are  capable  in  predicting  the  linear  and  nonlinear 
responses. 


Figure  35  Axial  creep  strain  (833)  and  transverse  ereep  strain  (s  n)  of  PZT-7A/LaRC-Si  eomposites  with 
19.6%  fiber  volume  eontents  due  to  a  eonstant  stress  033=10  MPa 
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Electric  flux  (D3  lO'^  C/m^) 


Figure  36  Electric  flux  along  the  polarization  axis  (D3)  of  PZT-7A/LaRC-Si  composites  with  19.6%  fiber 

volume  contents  due  to  a  constant  stress  a  33=10  MPa 
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Figure  37  Field  variables  of  the  aetive  PMC  panel  during  transient  heat  eonduetion  proeess 


41 


Figure  38  Comparison  of  electro-mechanical  properties  of  various  micromechanics  models  with 
experimental  data  (Chan  and  Unsworth,  1989)  for  the  PZT-7A/Araldite  D  fibrous  composite  with  fiber 
volume  fraction. 


Figure  39  Comparison  of  piezoelectric  constant  of  various  micromechanics  models  with  experimental 
data  (Furukawa  et  ah,  1976)  for  the  PZT-5/epoxy  particulate  composite. 
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Figure  40  Comparisons  of  effective  responses  from  MT  and  UC  models  for  active  fiber  composite  with 
fiber  volume  fraction  0.5  under  and  applied  electric  field  in  the  fiber  direction 


b)  Nonlinear  polarization  switching  responses 

The  fiber-unit-cell  and  particle-unit-cell  models  are  formulated  to  determine  the  effective 
hysteretic  polarization  and  butterfly  strain  responses  of  1-3  and  0-3  ferroelectric  composites, 
respectively.  A  phenomenological  constitutive  model  that  accounts  for  the  rate-dependent 
polarization  under  coupled  mechanical  and  electrical  effects,  proposed  by  Sohrabi  and  Muliana 
(2013)  and  discussed  in  Section  4.2.2,  is  used  for  the  ferroelectric  constituent  while  the  matrix  is 
considered  as  a  linear  elastic  or  viscoelastic  solid.  An  incremental  formulation  is  used  for  the 
constitutive  relations  and  unit-cell  models.  In  order  to  account  for  history-dependent  response 
from  the  time-dependent  ferroelectric  and  viscoelastic  matrix,  history  variables  are  introduced. 
The  linearized  constitutive  relation  for  each  constituent  is  written  as: 
a!  =0'AT' 


where  the  superscript  t  -  At  denotes  previous  time  and  S'  is  the  history  variables  of  the 
dependent  field  variables  S' .  The  components  of  O'  and  S'“'^'  for  the  ferroelectric  ceramics  and 
viscoelastic  materials  can  be  found  in  Lin  and  Muliana  (2014). 


Using  the  volume  averaging  scheme,  the  linearized  effective  field  variables,  denoted  by  an 
overbar,  of  the  UCs  at  current  time  t  are: 


S'  =  '^V^'"^E 


(m),t 


m=\ 


(24) 


where  N  is  the  number  of  subcells  (e.g.  AM-  for  a  fiber  UC  and  N=%  for  a  particle  UC).  The 

variable  S*'"^  '  is  the  average  field  variable  within  each  subcell  with  a  volume  fraction  . 
Following  a  linearized  constitutive  relation  for  each  constituent  in  Eq.  (23),  a  linearized 
constitutive  relation  for  the  overall  response  of  composites  at  current  time  Us  written  as: 

S' =0'AT'+S'"^'  (25) 
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In  order  to  relate  the  effective  incremental  independent  field  variable  in  the  UC,  AT' ,  to  the 
corresponding  incremental  field  variable  in  its  subcells,  ' ,  a  concentration  matrix  '  and 
a  vector  of  history  variables  '  at  current  time  t  are  defined  as  follows: 

^  ^  ^{my  _  ^26) 


Using  Eq.  (26)  and  (24),  the  effective  field  variables  in  Eq.  (24)  are  written  as: 

_  N  _  ^ 

m=\  m=\ 

_  N  _  ^ 

Q?  _  ^  ^t-At  _  ^  y(m)  ^{m\t-At  ^ 


(27) 


m=\ 


m=\ 


In  order  to  evaluate  the  concentration  matrices  and  history  variables  '  and  the 
constitutive  relations  for  the  different  constituents  together  with  the  linearized  micromechanical 
relations  from  the  fiber  and  particle  model  are  used  (see  Lin  and  Muliana,  2014). 


The  predictions  of  the  effective  hysteretic  polarization  and  butterfly  strain  responses  of  active 
fiber  composites  are  compared  with  existing  experimental  data  reported  by  Jayendiran  and 
Arockiarajan  (2013).  PZT-5Al/epoxy  fiber  composites  having  various  PZT-5A1  fiber  volume 
fractions  (VFs)  were  studied.  The  composites  were  subjected  to  cyclic  electric  field  with 
amplitude  of  =  ±2  MV/m,  along  X3  direction.  The  material  properties  of  the  PZT-5A1  are 

first  calibrated  from  the  experimental  data  of  PZT-5A1  (100%  fiber  content)  under  a  cyclic 
electric  field,  shown  in  Figures  41a  and  41b.  Jayendiran  and  Arockiarajan  (2013)  also  reported 
the  experimental  results  for  the  pure  PZT-5A1  subject  to  both  a  cyclic  electric  field  and  a 
constant  mechanical  stress  (i.e.,  0^33  =  -15  MPa)  shown  in  Figure  42  (circles).  The  existence  of 
the  compressive  stresses  influences  the  polarization  response  of  materials.  The  coercive  electric 
field  under  a  constant  compressive  stress  15  MPa,  =0.76  MV/m,  is  determined  when  the 

polarization  is  zero  as  shown  in  Figure  41a.  From  the  experimental  evidences  the  coercive 
electric  fields  vary  with  the  compressive  stresses.  Detailed  material  parameter  characterizations 
are  given  in  Lin  and  Muliana  (2014).  Micromechanical  model  predictions  of  composite  with  VF 
=  0.8  subjected  to  both  a  cyclic  electric  field  and  various  constant  mechanical  stresses  are  shown 
in  Fig.  43.  Figure  44  presents  responses  under  a  constant  compressive  stress  and  cyclic  electric 
fields  for  composites  with  different  fiber  contents.  The  agreements  on  the  effective  hysteretic 
polarization  between  the  micromechanical  predictions  and  experimental  data  are  good  for  each 
compressive  stress  while  the  estimations  of  the  effective  butterfly  strain  responses  show  some 
discrepancies  with  the  experimental  results.  The  height  of  the  butterfly  strain  response  of  the 
experimental  measurements  of  the  active  composite  is  larger  than  the  one  of  the  pure  PZT-5A1. 
One  would  expect  that  the  height  of  the  hysteretic  polarization  (saturated  polarization)  decreases 
as  the  volume  content  of  the  active  fiber  decreases  and  similarly  for  the  height  of  the  strain  in  the 
butterfly  strain  curves. 
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Figure  41  a)  Hysteretic  polarization  and  b)  butterfly  strain  responses  for  the  stress  free  PZT-5A1  due  to  a 
eyelie  eleetrie  field. 


Figure  42  a)  Hysteretie  polarization  and  b)  butterfly  strain  responses  for  the  PZT-5A1  due  to  a  eyelie 
eleetrie  field  and  a  eonstant  eompressive  stress. 
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Figure  43  Comparison  of  micromechanical  predictions  (solid  and  dashed  lines)  to  experimental  data 
(circles)  of  Jayendiran  and  Arockiarajan  (2013)  for  the  hysteretic  polarization  and  butterfly  strain 
responses  for  the  PZT-5Al/epoxy  1-3  active  composite  with  PZT-5A1  fiber  VF  =  0.8  subjected  to  both  a 
cyclic  electric  field  and  various  constant  mechanical  stresses. 
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Figure  44  Comparison  of  micromechanical  predictions  (solid  and  dashed  lines)  to  experimental  data 
(circles)  of  Jayendiran  and  Arockiarajan  (2013)  for  the  hysteretic  polarization  and  butterfly  strain 
responses  for  the  PZT-5A1 /epoxy  1-3  active  composite  with  various  PZT-5A1  fiber  VFs  subjected  to 
both  a  cyclic  electric  field  and  a  constant  mechanical  stress. 
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Parametric  studies  are  conducted  in  order  to  understand  the  electric-field  transfers  in  the  particle 
composites  as  the  active  inclusions  are  not  necessarily  in  direct  contact  with  the  electrodes, 
unlike  in  active  fiber  composites.  Particle  composites  with  epoxy  and  silver  matrix  are 
considered.  Epoxy  has  relatively  low  dielectric  constants,  while  silver  is  a  conductive  material. 
Figures  45  and  46  illustrate  the  polarization  and  strain  responses  in  particle  composites  with 
epoxy  and  silver  matrix,  respectively,  under  cyclic  electric  fields  and  a  constant  compressive 
stress.  As  seen  in  Fig.  45,  very  small  values  of  effective  polarization  and  strain  are  generated  in 
the  composites,  which  are  due  to  a  low  dielectric  matrix  (non-conductive),  and  as  a  result  the 
epoxy  matrix  is  unable  to  pass  significant  amount  of  charges  to  the  active  ferroelectric  particles. 
When  silver  (conductive)  matrix  is  used,  significant  hysteretic  polarization  and  butterfly  strain 
responses  are  observed  in  the  composites,  even  when  high  compressive  stress  that  can  suppress 
the  polarization  is  prescribed  to  the  composite.  This  shows  that  high  dielectric  constant  of  the 
conductive  matrix  causes  the  electric  charges  easily  pass  through  the  conductive  matrix  and 
eventually  large  amount  of  electric  fields  reaches  the  PZT-5A1  particles. 
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Figure  45  The  effeetive  polarization  and  strain  responses  for  the  stress  free  PZT-5A1 /epoxy  0-3  aetive 
eomposite  with  matrix  dieleetrie  eonstant  of  0.04x10“^  F/m. 
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Figure  46  Effective  hysteretic  polarization  and  butterfly  strain  responses  for  the  PZT-5A1 /silver  0-3 
active  composite  with  PZT-5A1  particle  undergoing  both  a  cyclic  electric  field  and  a  constant  mechanical 
stress  ^33  =  -45  MPa. 


c)  Response  of  active  functionally  graded  materials 

A  piezoelectric  bimorph  beam,  which  is  formed  by  bonding  two  thin  layers  of  piezoelectric 
ceramic  plates,  is  often  used  as  actuators.  An  actuating  mechanism  is  induced  by  applying 
electric  fields  through  the  thickness  of  the  piezoceramics  plate  so  that  one  of  the  piezoelectric 
ceramic  plates  would  elongate  in  the  longitudinal  direction  and  the  other  plate  would  contract  in 
the  longitudinal  direction.  The  relatively  large  differences  in  the  deformations,  i.e.,  tension  and 
compression,  in  the  longitudinal  direction  of  the  two  piezoelectric  ceramic  plates  lead  to  high 
stress  discontinuities  at  the  interfaces  between  the  piezoelectric  plates.  The  high  stress 
discontinuities  at  the  interface  can  cause  debonding  and  eventually  failure  of  the  piezoelectric 
bimorph  beam.  Functionally  graded  piezoelectric  beams  (FGPBs)  have  been  considered  as  one 
of  the  promising  solutions  to  relieve  the  high  stress  discontinuities  at  the  interface  of  the 
piezoelectric  ceramic  plates  by  introducing  gradual  changes  in  the  compositions  of  the 
piezoelectric  ceramic  plates,  see  for  example,  Wu  et  al.  (1996),  Jin  and  Meng  (2003),  Takagi  et 
al.  (2002). 
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This  study  presents  analyses  of  the  nonlinear  electro-mechanical  response  of  polarized  FGPBs 
and  polarization  switching  behavior  in  FGPBs,  undergoing  small  deformation  gradients.  The 
lateral  deformations  in  the  FGPBs  are  due  to  applications  of  electric  fields  through  the  thickness 
of  the  beams.  The  functionally  graded  beam  is  discretized  into  several  graded  layers  through  its 
thickness  as  shown  in  Fig.  47.  The  interfaces  between  layers  are  assumed  perfectly  bonded  so 
that  the  charge  continuously  flows  through  the  thickness,  i.e.,  charges  in  every  layer  of  the 
FGPBs  are  assumed  uniform.  Each  layer  is  comprised  of  different  compositions  of  the  active 
(piezoelectric)  inclusions  and  conductive  matrix,  i.e.,  metallic  matrix.  The  microstructure  in  each 
layer  is  idealized  as  solid  spherical  inclusions  randomly  distributed  in  a  homogeneous  matrix  and 
it  is  assumed  that  all  particles  are  fully  surrounded  by  homogeneous  matrix  without  any  contents 
of  pores  and  defects.  A  unit-cell  model  consisting  of  eight  inclusion  and  matrix  subcells  is  then 
considered  for  obtaining  the  overall  (homogenized)  electro-mechanical  properties  and  responses 
in  each  layer.  The  inclusions  and  the  matrix  are  assumed  perfectly  bonded  and  the  spatial 
variations  of  the  field  variables  in  each  subcell  are  assumed  uniform.  Further  discussion  on  the 
active  functionally  graded  beams  can  be  found  in  Lin  and  Muliana  (2014b). 
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Figure  47  Schematic  of  a  functionally  graded  piezoelectric  beam  and  its  microstructural  approximation. 


The  predictions  of  the  tip  deflections  and  curvatures  of  the  piezoelectric  beams  are  first 
compared  with  existing  experimental  data,  which  is  limited  to  responses  of  polarized 
piezoceramics.  Parametric  studies  are  also  performed  to  study  the  effects  of  thickness  of  the 
beam,  constituent  volume  fraction,  and  excitation  frequency  on  the  overall  performance  of  the 
beams  under  cyclic  electric  field  inputs.  First,  experimental  data  reported  by  Alexander  and  Brei 
(2005)  are  used  for  comparisons.  The  FGPB  consists  of  PZT-856  and  BT  constituents  with  five 
piezocomposite  layers,  i.e.,  pure  PZT-856,  0.95/0.05,  0.90/0.10,  0.85/0.15,  and  0.80/0.20  VF  of 
PZT-856/BT.  Two  sides  of  the  beams  were  electroded  with  1500  A  of  AuPd.  Two  possible 
microstructures  of  the  PZT-856/  barium  titanate  functionally  graded  beams  are  considered  since 
their  detailed  microstructural  morphologies  are  not  discussed  in  the  literature.  The  first 
microstructure  considers  PZT-856  inclusions  dispersed  in  homogeneous  barium  titanate  matrix, 
in  which  the  PZT  inclusions  are  fully  surrounded  by  barium  titanate.  The  second  microstructure 
is  the  inverse  of  the  first  one.  Figure  48  shows  the  comparisons  between  the  predictions  and 
experimental  measurements  of  the  tip  deflection  of  the  studied  functionally  graded  beams.  The 
predictions  from  the  PZT-856  inclusions  dispersed  in  the  barium  titanate  matrix  (solid  lines)  are 
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in  good  correlations  with  experimental  results.  The  predictions  of  the  piezocomposite  with 
barium  titanate  inclusions  dispersed  in  the  PZT-856  matrix  (dashed  lines)  underestimate  the  tip 
deflection.  Another  functionally  graded  piezoelectric  beam  considered  here  was  developed  by 
Takagi  et  al.  (2002),  which  consists  of  platinum  (Pt)  particles  dispersed  in  polarized  PZT 
matrices.  Figure  49  depicts  the  predictions  and  experimental  data  of  the  curvature  of  the 
functionally  graded  piezoelectric  beam  against  the  applied  voltage.  It  is  expected  that  for  the  low 
voltage  linear  theory  works  well. 


-60' - ^ ^ ^ - ' 
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Figure  48  Relationships  between  a  total  applied  voltage  and  the  tip  deflection  of  the  five-layer 
functionally  graded  piezoelectric  beam. 


Figure  49  Relationships  between  a  total  applied  voltage  and  the  curvature  of  the  functionally  graded 
piezoelectric  bender.  The  solid  (Nonlinear  PZT)  and  dashed  (Linear  PZT)  lines  show  the  predictions  from 
the  model  while  solid  circle  symbols  are  experimental  measurements  of  curvatures. 
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Parametric  studies  are  now  conducted  for  a  functionally  graded  piezoelectric  cantilever  beam 
with  five  piezocomposite  layers  undergoing  cyclic  electric  fields  in  order  to  examine  the  effect  of 
polarization  switching  response  on  the  tip  deflections.  Each  layer  has  different  compositions  of 
PZT-51  particle  and  Platinum  matrix.  Higher  frequency  results  in  smaller  hysteretic  responses,  as 
seen  in  Fig.  50.  This  is  because  the  PZT-51  constituent  does  not  have  enough  time  to  experience 
“creep-like”  behavior  at  higher  frequency. 


Voltage  (V) 

Figure  50  Effect  of  the  excitation  frequency,  /  =  a)  0.5,  b)  1,  and  c)  10  Hz,  on  the  tip  deflection 


4.3.3)  Micromechanical  models  with  detailed  microstructural  morphologies 

Micromechanics  models  with  idealized  and  simplified  shapes  of  inhomogeneities,  such  as  the 
ones  discussed  in  Sections  4.3.1  and  4.3.2,  have  been  widely  used  to  obtain  the  average 
(macroscopic)  mechanical  response  of  different  composite  materials.  One  of  the  purpose  of  this 
study  is  to  examine  whether  the  composites  with  irregular  shapes  of  inhomogeneities,  such  as  in 
the  aluminum-alumina  (AI-AI2O3)  and  barium  titanate-silver  (BaTiOs-Ag)  composites,  can  be 
approximated  by  considering  idealized  and  simplified  shapes  of  inhomogeneities  in  determining 
their  overall  macroscopic  mechanical  responses.  It  is  noted  that  AI-AI2O3  composites  are  often 
used  as  host  structures  in  high  temperature  application,  while  BaTiOs-Ag  composites  are 
promising  for  active  composites  used  at  elevated  temperatures. 
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For  this  purpose,  two  composite  systems,  i.e.,  AI-AI2O3  and  BaTiOs-Ag,  are  manufactured  using 
powder  metallurgy  method  and  characterized  for  their  thermo-mechanical  and  thermo-electro- 
mechanical  properties. 

a)  AI-AI2O3  composite  system 

AI-AI2O3  composites  with  different  AI2O3  contents  and  particle  sizes  have  been  manufactured 
and  studied.  Information  on  processing  of  the  AI-AI2O3  composites  can  be  found  in  Gudlur 
(2013).  The  effects  of  shapes  of  the  inhomogeneities  (inclusions  and  pores)  that  are  randomly 
dispersed  in  a  homogeneous  ductile  matrix  on  the  overall  thermo-mechanical  response  of  Al- 
AI2O3  composites  are  examined.  The  micromechanical  models  with  detailed  microstructural 
morphologies  are  constructed  from  the  microstructural  images  of  the  composites  obtained  from 
the  scanning  electron  microscope  (SEM).  The  SEM  images  of  the  composites  are  converted  into 
finite  element  (FE)  meshes  using  software  OOF2*  and  ABAQUS  FE  code  is  used  to  analyze  the 
overall  response  of  the  AI-AI2O3  composite.  Figure  51  shows  the  SEM  image  of  a  composite 
with  20vol%  AI2O3,  where  the  alumina  particles  appear  light  in  color,  aluminum  matrix  in  grey, 
and  voids/pores  appear  dark  in  color.  The  size  of  the  microstructural  image  is  217pm  x  173  pm. 
This  micrograph  image  is  divided  into  uniform  sub-images  (regions).  As  an  example,  the  image 
is  divided  into  12  uniform  regions  with  50pm  x  50pm.  Four  different  50  pm  x  50  pm  square 
regions  having  different  microstructures  are  randomly  chosen  as  representative  microstructures 
of  the  composite.  Four  representative  volume  elements  (RVEs)  which  are  referred  as  FE  meshes 
1,  2,  3,  and  4  are  randomly  selected.  The  overall  linear  elastic  modulus  of  the  composites  from 
the  selected  RVEs  is  first  studied.  Convergence  studies  in  terms  of  number  of  elements  used  in 
FE  analyses  and  also  the  sizes  of  RVE,  in  which  the  RVE  sizes  are  varied  by  selecting  two 
regions  with  100pm  x  100pm,  are  performed.  The  two  RVEs  are  referred  as  FE  meshes  5  and  6. 
The  elastic  moduli,  determined  from  different  RVE  sizes  and  number  of  elements,  are  compared. 
The  FE  simulation  is  carried  out  in  two  steps.  The  first  step  is  prescribing  a  uniform  temperature 
change  from  200°C,  which  is  considered  as  the  stress-free  temperature  based  on  a  previous  work, 
Gudlur  et  al.  (2014),  to  25°C  and  examining  the  residual  (remaining  thermal)  stresses. 

The  linear  responses  from  the  FE  analyses  are  first  studied  and  compared  to  the  ones  obtained 
from  the  experiment  using  nondestructive  testing,  resonant  ultrasound  spectroscopy  (RUS). 

Table  10  presents  the  values  of  elastic  moduli  and  Poisson’s  ratios  from  the  six  FE  meshes 
considered.  The  first  four  FE  meshes  (FE-1,  FE-2,  FE-3,  and  FE-4)  have  50pm  x  50pm  square 
region,  while  FE-5  and  FE-5  have  100pm  x  100pm  square  region.  The  average  size  of  the 
element  in  the  FE  analyses  with  coarse  elements  is  3.27x10'"^  x  3.27x10'"^  mm^,  while  the  one 
with  finer  elements  is  l.OhxlO  '^x  l.OhxlO  '^mm^.  The  coarse  elements  are  capable  in  capturing 
the  overall  properties  of  the  composites  while  significantly  reduce  the  computational  cost.  It  is 
also  seen  that  the  smaller  RVE  size  is  capable  in  representing  the  heterogeneities  in  composites. 
For  comparison,  the  overall  elastic  properties  of  composites  without  considering  the  thermal 
(residual)  stress  effect  are  shown  in  Table  11.  The  properties  from  four  FE  meshes  are  compared 
to  the  Mori-Tanaka  micromechanics  model.  It  is  seen  that  the  elastic  moduli  without  the  residual 
stress  effect  are  about  12  GPa  higher  than  the  ones  with  the  residual  stress  effect.  It  is  also  noted 
that  the  effective  elastic  properties  from  the  four  FE  analyses  are  closed  to  the  ones  of  Mori- 
Tanaka  model  with  20%  alumina  volume  contents,  in  which  an  idealized  microstructure  of  a 
composite  with  a  solid  sphere  embedded  in  homogeneous  matrix  is  considered.  From  the  above 
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analyses  one  ean  observe  that  the  existence  of  residual  (thermal)  stresses  in  the  constituents 
significantly  influence  the  macroscopic  elastic  properties  of  the  composites. 


FE  models 


Figure  5 1  a)  SEM  Image  of  AI-AI2O3  composite  sample  with  20%  alumina  contents;  b)  Four  random  square 
micrographs  of  the  composite  sample  and  their  corresponding  FE  microstmctural  models  1,  2,  3,  and  4 
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Table  10.  Area  Fraetion  of  Alumina  (%AF),  Porosity  (%Pores),  elastie  modulus,  and  Poisson’s  ratio  of 
the  FE  models  generated  from  eomposite  with  20%  Volume  Fraetion  (%VF)  of  alumina.  For  eomparison, 
results  of  the  experiment  from  the  nondestruetive  testing,  resonant  ultrasound  speetroseopy  (RUS),  are 
also  provided  in  the  table. 


Model 

%AF  of  Alumina 

%Pores 

E  (GPa) 

V 

FE-1 

20.2 

0.4 

71.23770.82'^ 

0.397/0/402 

FE-2 

21.7 

0.1 

72.57/72.15 

0.393/0.399 

FE-3 

20.9 

0.7 

69.44/69.06 

0.390/0.395 

FE-4 

21.2 

0.0 

75.53/75.08 

0.393/0.395 

FE-5 

21.7 

0.3 

71.43 

0.405 

FE-6 

19.1 

0.8 

70.03 

0.396 

RUS 

22.6 

0.5 

75.15 

0.310 

"using  coarse  meshes;  ’’using  finer  meshes 

Table  1 1  Effeetive  elastie  properties  of  composites  with  20%  VF  of  alumina,  without  residual  (thermal) 
stress  effect 


Model 

E  (GPa) 

V 

FE-1 

83.28 

0.313 

FE-2 

86.67 

0.311 

FE-3 

83.70 

0.314 

FE-4 

88.05 

0.315 

Mori  Tanaka 

87.00 

0.330 

The  second  micromechanics  model  is  constructed  by  considering  idealized  shapes  of  the 
inhomogeneities,  while  maintaining  the  same  volume  contents  and  locations  of  the 
inhomogeneities  as  the  ones  in  the  micromechanics  model  with  microstructural  details.  In  this 
case,  the  overall  inelastic  and  degradation  responses,  incorporating  the  thermal  stress  effect,  from 
the  detailed  microstructural  characteristics  is  compared  to  the  ones  determined  with  idealized 
microstructures.  The  alumina  is  modeled  as  linear  elastic  while  the  aluminum  is  assumed  to 
exhibit  an  elastic-plastic  deformation  following  the  Hencky-Ilyushin  model.  In  order  to 
incorporate  the  degradation  in  the  ceramics  caused  by  relatively  high  localized  tensile 
stresses/strains,  the  following  rate  degradation  model  for  the  ceramics  is  used  j  =  (i- 

where /is  the  general  function  of  strains.  The  elastic  property,  i.e.,  elastic  modulus,  of  the 
ceramics  changes  with  the  degradation,  E{d).  It  is  noted  that  this  study  does  not  aim  at  rigorously 
model  the  degradation  in  the  composite  constituents,  instead  it  is  done  to  briefly  examine  the 
overall  macroscopic  response  of  composites  in  case  of  some  localized  degradation  occurs  in  the 
constituents.  Figure  52  illustrates  the  stress  and  degradation  fields  in  FE  meshed  model  #1.  As 
seen  in  the  figures,  the  composite  with  actual  shapes  of  inhomogeneities  shows  stress  and 
degradation  fields  with  higher  magnitude  as  compared  to  the  composite  with  circular  shapes  of 
inhomogeneities,  especially  in  the  brittle  inclusions.  It  is  noted  that  sharp  edges  of  the 
inhomogeneities  in  the  microstructures  with  actual  shapes  of  inhomogeneities  lead  to  high 
localized  stresses.  Figure  53  presents  the  corresponding  overall  responses  from  the  two  micro¬ 
mechanics  models,  which  show  very  close  stress-strain  behaviors. 
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Figure  52  Transverse  stress  and  degradation  eontours  from  FE  meshed  model#!  with  aetual  (top)  and 
idealized  (bottom)  inhomogeneities 
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Figure  53  Maeroseopie  (effeetive)  stress-strain  responses  of  the  AI-AI2O3  eomposite  with  aetual  and 
idealized  mierostruetures. 
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b)  BaTiOs-Ag  composite  system 

BaTiOs-Ag  composites  with  different  silver  contents  have  been  manufactured  and  studied.  Some 
information  on  the  manufacturing  and  testing  of  BaTiOs-Ag  composites  can  be  found  in  Xing  et 
al.  (2013).  Similar  to  the  AI-AI2O3  composites,  micromechanics  models  of  BaTiOs-Ag 
composites  with  microstructural  details  are  formulated,  see  Fig.  54  as  an  example.  The 
micromechanics  models  are  used  to  study  the  thermal,  electrical,  and  mechanical  properties  of 
the  composites  at  different  temperatures:  -50°C  to  250°C.  Phase  transformations  in  the  BaTiOs 
are  observed  around  25°C  and  125°C.  Figure  54  shows  the  thermal  expansion  coefficient  of 
BaTiOs-Ag  composites  with  different  silver  contents  and  at  three  different  phases,  while  Fig.  56 
illustrates  the  thermal  conductivity  at  room  temperature.  The  micromechanical  models  are 
capable  in  capturing  the  overall  thermal  properties  of  composites.  Similar  predictions  are  done 
on  the  mechanical  and  electrical  properties. 


Figure  54  Example  of  microstmctural  model  of  BaTiOs-Ag  with  10%  silver  content 
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Figure  55  Thermal  expansion  eoeffieient  of  BaTiOs-Ag  at  different  phases  and  silver  eontents 
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Figure  56  Thermal  eonduetivity  response  of  BaTiOs-Ag  at  room  temperature 


4.4)  Structural  Analyses 

Multi-level  homogenization  schemes  that  incorporate  different  constitutive  material  models  of 
the  constituents  and  microstructural  details  have  been  formulated  for  predicting  the  effective 
response  of  structures  comprising  of  heterogeneous  materials.  In  multi-scale  models,  numerical 
methods  are  often  used  to  find  solutions  to  boundary  value  problems  (BVPs)  at  multiple  length- 
scales.  In  this  project,  a  multi-scale  model  is  formulated  based  on  an  integrated  simplified 
micromechanical  model  and  FE  framework,  which  is  used  to  analyze  nonlinear  time-dependent 


thermo-electro-mechanical  responses  of  smart  structures.  The  smart  structures  can  comprise  of  a) 
composite  host  structures  integrated  with  active  materials  made  of  piezoceramics  or  active 
composites,  b)  arrangements  of  active  and  inactive  components  of  various  materials.  The 
purpose  is  to  perform  shape  changes  in  the  smart  structures  controlled  by  the  electro-active 
components. 

In  this  study,  coupled  transient  heat  conduction  and  thermo-electro-mechanical  responses  are 
considered.  Finite  element  (FE)  method  is  used  to  obtain  solutions  to  BVPs  in  the  active 
structures.  It  is  assumed  that  the  heat  conduction  influences  the  electro-elastic  response,  but  the 
electro-elastic  deformations  do  not  influence  the  heat  conduction  in  the  active  materials, 
allowing  for  sequentially  coupled  heat  conduction  and  deformation  analyses.  When  fully  coupled 
heat  conduction  is  considered,  user  subroutine  materials  UMAT  and  UMATH  within  ABAQUS 
FE  are  used  in  order  to  determine  the  solutions  to  field  variables,  see  Khan  and  Muliana  (2012b). 


The  multi-scale  analyses  of  heat  conduction  and  thermo-electro-mechanical  deformations  within 
3D  continuum  finite  elements  are  briefly  summarized  as  follows.  Let  It”1  =  |r',r^,...,r^'^lbe 


the  nodal  temperatures  in  a  continuum  element  with  number  of  nodes  Nd.  An  overbar  indicates 
the  macro-scale  quantities  and  T  denotes  transpose  of  a  vector  or  matrix.  The  temperature 
gradients  at  the  macro-level,  which  are  sampled  at  material  (Gaussian)  points,  are  given  as: 


dT 

dXj 


or  (p  =  B^T 


(28) 


where  is  the  gradient  matrix.  Using  the  principle  of  virtual  work  the  energy  balance  equation 
when  only  heat  conduction  is  considered,  is: 

j  n’^CvNJFT"  +  J  =  J  +  J  ^'^JdA 

V  V  V  A 

(29) 

C^t"+K^T"  =R'‘+R^; 

dt 

where  N  is  the  shape  function  matrix,  Cv  is  the  macro-scale  heat  capacity;  ky  is  the  component 

of  the  macro-scale  (effective)  thermal  conductivity  ;  and  r  and  /  are  the  rate  of  internal  heat 
generation  per  unit  volume  and  surface  heat  flux,  respectively;  V  is  the  volume  of  the  body,  and 
A  is  the  surface  area  where  the  heat  flux  is  prescribed.  A  general  element  should  be  able  of 
incorporating  any  orientation  in  a  global  coordinate  system.  This  can  be  done  through  rotational 
tensor  transformations.  Macro-scale  thermal  and  physical  properties,  i.e.,  ky  and  Cv  are  needed 


— T  — T 

in  order  to  form  C  and  K  matrices.  When  homogeneous  materials  are  considered,  the  heat 
conductivity  and  heat  capacity  of  the  homogeneous  materials  can  be  used  directly.  In  case 
composite  materials  are  considered,  the  effective  properties  discussed  in  Section  4.3.1  are  used. 
It  is  noted  that  if  heat  conduction  through  an  active  body  causes  significant  temperature  changes 
in  the  body,  variations  of  the  thermal  conductivity  and  heat  capacity  with  temperatures  are  not 
negligible,  leading  to  a  nonlinear  problem.  In  order  to  obtain  solutions  for  the  temperature  field 
at  the  macro-scale,  backward  difference  method  for  the  temporal  integration  and  Gaussian 
quadrature  method  for  the  spatial  integration  are  adopted.  At  each  time  increment,  a  linearized 
relation  in  Eq.  (29)  is  used  as  a  starting  point  for  obtaining  the  temperature  field  at  current  time 
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and  the  Newton-Raphson  iterative  method  is  then  used  to  eorrect  for  errors  due  to  linearization. 
Let  At  be  an  incremental  time  and  the  temperature  time  derivative  at  current  time  (t  +  At)  is 

approximated  as:  TP+At  ~  - — ;  where  T”  is  the  nodal  temperature  at  previous  time  t  and 


-.n 

0 


refers  to  initial  conditions  for  the  nodal  temperatures.  If  the  internal  heat  generation  is  absent 


and  the  surface  heat  flux  is  the  external  stimuli,  then  .  With  material  properties  vary 

with  temperatures,  the  linearized  equilibrium  equation  for  an  entire  structure  is  now  written  as: 


^  T  rp  I  T  T'  ^ 

^  t+At  t+At  t+At  t+At  ~  ^  t+At 


(30) 


A  multi-scale  time-integration  algorithm  is  formulated  to  obtain  solutions  for  nodal 
temperatures  in  Eq.  (28),  which  is  summarized  as  follows: 

1)  Input  variables:  At ,  ,  K[ ,  Cf ,  and 


2)  Calculate 


(linearized  solution,  iteration  counter  A=0): 


1 

— C 

\At 


rTiW,0 

^t+At 


_  u  ^  T  rji  n 

~  ^t+At 

At 


3)  Iterate  for  A=0,I,2,  . . . 

3.1)  Evaluate:  7)+^^  =  where  (m)  indicates  subcell’s  number; 

t  ,  t(..)  );  from  Eq.  (3.5) 

C„f,^  from  Eq.  (3.6) 
C,’;‘.d,“‘iK"3,fi'»”Eq.(3.2) 


3.2)  Define  residual  vector: 
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In  order  to  define  the  above  Jaeobian  matrix,  the  variations  of  the  maero-seale  thermal 
eonduetivity  and  heat  eapaeity  with  nodal  temperatures  are  needed: 


^^t+At  _  yi^) 
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4)  Output  variables:  T^"  ,  R^t+At  ^  ^t+At 
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The  above  time-integration  algorithm  is  formulated  to  solve  for  the  nonlinear  heat  eonduetion 
equations  at  multiple  length  seales:  homogeneous  eonstituents,  mieroseopie  structures  of  active 
composites,  and  global  active  structures. 


Similar  to  the  heat  conduction  analyses,  the  corresponding  thermo-electro-elastic  response  of  the 
structure  subjected  to  external  mechanical  and  non-mechanical  stimuli  can  be  determined.  Let 


(—1  —1  —1  —Nd  —Nd  —Nd 
—  I  ^1  ^  ^  ^3  5  """5  U-^  ?  ^2  5  ^3 


and 


(—1  —2  —Nd 

=  W  ,(p 


be  the  nodal  displacement  and 


electric  potential  vectors,  respectively.  The  strain  and  electric  field  at  the  macro-level  are  defined 
as: 


^=BrjntUl 

or  s  =  B"u” 

(31-1) 

dxt 

or  E  =  B^o" 

(31-2) 

where  B“  and  B***  are  the  spatial  derivative  of  shape  functions  related  to  the  macro-scale  strain 
and  electric  field,  respectively.  The  energy  balance  equations  related  to  the  mechanical  and 
electrical  components  in  linearized  forms  are  (in  absences  of  body  forces  and  body  charges): 


(32-1) 


j  B"^  CB“u"-e^B^a>"-pNT"  = 


j(B^'^  ^“U"+i®^®"+pNT" 
+K^O”  =-F^^  -l-R^ 


JjF  =  j 


(32-2) 


where  t(  and  q  ^  denote  the  components  of  the  surface  traction  and  surface  charge,  respectively, 
prescribed  at  the  macro-scale.  The  components  of  the  electro-elastic  coupling  and  dielectric 

stiffness  matrices  depend  on  the  potential  gradient  at  current  time.  The  nodal  temperatures 

T”  are  obtained  by  first  solving  the  heat  conduction  equation.  Since  this  study  consider 

nonlinear  time-dependent  thermo-electro-mechanical  responses,  the  linearized  relations  in  Eq. 
(32)  are  obtained  at  each  time  increment  in  order  to  provide  trial  solutions  of  the  nodal 
displacements  and  electric  potential  and  the  Newton-Raphson  iterative  method  is  then  used  to 
correct  for  errors  from  the  linearization.  The  multi-scale  analyses  at  current  time  ( t  +  At ) 
after  assembling  of  elements  are  summarized  as  follows. 

1)  Input  variables  are  At ,  ,  Kf“ ,  Kf ,  and 


2)  Evaluate  =  NT"  =  T, 
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V.“U  and  kJ 

3)  Calculate  and  for  k=Q  (trial  solutions): 
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Struetural  analyses  are  performed  in  order  to  simulate  the  overall  performance  of  several  active 
structures  subjected  to  various  histories  of  thermo-eleetro-meehanical  behaviors.  Analytical  and 
numerical  methods  have  been  used. 


a)  Shape  changes  in  active  beams  due  to  electric  field  inputs 

Deformations  in  a  fiber  reinforced  laminated  eomposite  plate  bonded  with  four  piezoelectrie 
patches  due  to  eleetric  field  inputs  are  analyzed.  The  eomposite  plate  is  first  assumed  linear 
elastie  with  the  following  elastie  material  properties:  £'=36000  MPa,  v=0.25  in  the  longitudinal 
fiber  axis.  The  plate  is  clamped  along  one  of  its  50  mm  side  surfaee  and  the  PZT  patches  are 
bonded  perfectly  to  the  host  structure.  The  potential  on  the  surfaees  of  the  PZT  patehes  that  are 
in  eontaet  to  the  host  structure  is  grounded  to  zero.  The  PZT  patehes  are  uniformly  subjeeted  to  a 
potential  gradient  through  their  thickness  E^{t)  =  1.2sin(2;T  ft)  MV/m  so  that  they  experience 

expansion  in  their  planar  direetion,  thus  indueing  bending  to  the  eomposite  plate.  Figure  57 
depicts  the  lateral  displacement  measured  at  the  mid-section  of  the  free  end  due  to  the  input 
electrie  field  applied  uniformly  to  the  four  PZT  patehes,  which  show  vibration  of  the  eomposite 
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plate.  It  is  seen  that  the  highest  deformation  occurs  at  the  first  quarter  cycle,  when  the  electric 
field  reaches  1.2  MV/m  (point  A).  Upon  removal  of  the  electric  field,  the  remanent  deformation 
is  shown  (point  B)  and  when  the  electric  field  reaches  the  coercive  limit  of  the  PZT  patches  no 
deformation  is  shown  in  the  plate  (point  C)  due  to  depolarization  of  the  PZT  patches.  Continuing 
applying  the  electric  field  till  it  reaches  the  lowest  peak  -1.2  MV/m  (point  D)  trigger  bending  in 
the  composite  plate.  It  is  also  seen  that  due  to  the  time-dependent  effect,  the  highest 
displacement  is  at  the  first  quarter  cycle  and  after  several  cycles,  saturation  in  the  butterfly 
displacement  loop  is  achieved.  The  corresponding  deformed  shapes  of  the  composite  plate  at 
several  instant  of  times  during  the  first  cycle  are  illustrated  in  Fig.  58.  It  is  also  possible  to  create 
different  deformed  shapes  by  considering  various  histories  of  electric  field  inputs.  For  example, 
to  induce  twisting  of  the  plate  different  electric  field  inputs  can  be  applied  to  the  four  patches. 
The  tip-deflection,  illustrated  in  Fig.  59,  is  generated  by  prescribing  the  following  electric  field 
input  £'3(0  =  1.2sin(2;T  ft)  MV/m  to  the  two  PZT,  and  Eft)  =  1.2sin(2.;r  ft-n !  2)  MV/m  to  the 
other  two  patches. 

The  effect  of  time  dependent  response  of  a  viscoelastic  host  structures  with  piezoelectric  patches 
on  the  overall  mechanical  deflection  is  also  examined.  Figure  60  illustrate  the  lateral  deflection 
in  smart  beams  activated  by  applying  electric  field  to  the  patches.  The  patches  are  activated  with 
the  same  time  function  with  the  frequency  of  0.1  Hz.  Elastic  and  viscoelastic  host  structures  are 
considered.  It  is  seen  that  responses  of  the  lateral  deflections  of  the  elastic  and  viscoelastic  host 
structure  are  quite  different.  A  delayed  response  of  the  viscoelastic  host  structure  due  to  the 
electric  field  inputs  is  observed.  It  is  noted  that  small  loops  are  seen  in  the  response  at  the  highest 
amplitude  of  electric  field,  which  could  be  due  to  a  superimposed  of  two  different  hysteresis 
behaviors  from  the  mechanical  response  of  host  structure  and  cyclic  electric  field  of  the  patches. 
The  small  loops  near  the  peaks  of  the  hysteresis  responses  are  experimentally  observed  in 
piezoelectric  micro  fiber  composite  (MFC)  actuators  (Usher  et  al.,  2013)  and  are  also  observed 
as  distortion  of  hysteresis  loop  while  resistor  is  added  to  the  actuation  circuit  of  BaTiOs  (Jaffe 
2012). 


b)  Tip  displacement  Displacement 

EjiMV/m)  a)  Input  electric  field  (mm)  A 


Figure  57  Shape  ehanging  performanee  in  a  glass/fiber  laminated  eomposite  plate  subjeeted  to  a  hysteretie 
eleetrie  field  input.  The  deformation  is  eontrolled  through  the  use  of  piezoeleetrie  eeramies  patehes. 
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Time=0.5sec  (B) 


Time=0.66  sec  (D) 


0  I  0.075  10.15  mm 


Figure  58  Deformed  shapes  of  smart  beams  due  to  eleetrie  field  inputs  applied  to  the  piezoeleetrie 
patehes. 


Uj  tip  deflection  (mm) 
0.1? 


Figure  59  Tip  defleetion  of  the  eantilever  plate  due  to  a  non-uniform  eyelie  eleetrie  field  applied  to  the 
piezoeleetrie  patehes 
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Q  25  Lateral  dispi 


Figure  60  Lateral  defleetions  in  smart  eomposites  with  elastie  and  viseoelastie  host  struetures 

A  nonlinear  electro-mechanical  response  of  bimorph  beam  subjected  to  cyclic  electric  field  input 
with  frequency  of  IHz  is  also  considered,  following  an  experimental  test  reported  by  Li  and  Fang 
(2004).  The  beam  consists  of  three  layers  with  piezoelectric  layers  attached  to  the  top  and  bottom 
surfaces.  The  nonlinear  electro-mechanical  coupling,  discussed  in  Section  4.2.,  is  adopted,  which 
give  the  following  electro-mechanical  relation  (jjj  =  and  the  time- 

dependent  kernel  function  is  taken  as  K{t)  =  Kq+K^  exp”^' .  The  beam  is  under  simply 

supported.  Analytical  solution,  using  the  Laplace  transform  method,  is  considered  in  predicting 
the  bending  deformation  of  the  beam  due  to  the  electric  field  input.  Figure  61  illustrates  the 
electro-mechanical  response  of  the  beam  obtained  from  analytical  and  experimental  results. 

Good  correlations  between  the  model  and  experiment  are  observed. 


Figure  61  Analytical  and  experimental  results  of  bending  deformation  in  active  beam  due  to  electric  field. 


65 


Piezoelectric  materials  can  also  be  used  for  sensing  and  actuating,  such  as  detecting  and 
suppressing  the  deformation  in  a  beam.  Song  et  al.  (2004)  conducted  a  test  on  a  smart  beam  with 
PZT  sensor  and  actuator  (Fig.  62a),  in  which  the  PZT  is  used  to  sense  an  excessive  deformation 
in  the  beam  due  to  a  temperature  change  and  the  PZT  actuator  is  used  to  minimize  the  lateral 
deformation.  In  this  study,  the  suppression  of  the  lateral  deformation  in  the  beam  using  PZT 
actuator  is  simulated.  For  this  purpose,  the  placement  of  the  PZT  patches  is  slightly  modified 
(Fig.  62b).  The  steady  state  temperature  field  is  assumed  in  the  composite  beam.  Figure  63 
compares  the  numerically  evaluated  vertical  displacement  of  the  beam  due  to  temperature 
gradient,  for  both  the  actuated  and  the  unactuated  states.  The  tip  displacement  of  the  unactuated 
beam,  predicted  numerically,  is  1 1.46  mm  compared  to  1 1  mm  as  reported  in  Song  et  al.  (2004). 
Hence,  there  is  good  agreement  between  both  the  theoretical  predictions  and  experimental 
results.  The  numerically  calculated  actuating  voltage  is  45V  as  compared  to  150V 
experimentally.  The  difference  in  the  magnitudes  of  both  the  voltages  could  be  attributed  to  the 
modified  geometry,  since  the  length  of  actuator  in  the  modified  geometry  is  greater  than  that  of 
the  original. 


PZT  5A  PZT  5H 


2.79mm 


70*0.32  mm2 
46*0.254  mm2 


(a)  Composite  beam  with  .3  PZT  5H  actuators  and  one  PZT  5A  sensor 

PZT5H 


Rigid  clamp  849  mm 


■» 


2.79mm 


245*0.32  mm2 


(b)  Composite  beam  with  PZT  5H  surface  bonded  actuator 

Figure  62  A  smart  composite  beam  with  PZT  patches:  (a)  the  geometry  used  by  Song  et  al.  ;  (b)  the 
modified  geometry  used  in  the  current  numerical  study 
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Figure  63  Comparison  of  the  vertical  displacement  obtained  numerically  to  that  reported  in  Song  et  al. 
(2004)  for  the  beam  in  actuated  and  unactuated  states.  Curves  highlighted  by  (‘o’)  depict  the  experimental 

results 


Next,  the  integrated  micromechanical  and  FE  analyses  are  used  to  study  the  deformation  in  smart 
sandwich  composites  having  polyurethane  foam  core  and  piezoelectric  sensors  embedded  in  the 
laminated  facesheets  (Fig.  64).  The  simplified  micromechanical  model  for  fiber  composites 
discussed  in  Section  4.3  is  implemented  in  continuum  finite  element  and  used  to  obtain  the 
overall  response  of  facesheets  (see  Kim  et  al.  201 1  for  detailed  discussion).  Viscoelastic 
constitutive  models  are  considered  for  the  polyurethane  foam  and  epoxy  matrix  in  the  facesheet. 
Linear  electro-mechanical  response  is  assumed  for  the  piezoelectric  wafers.  It  is  noted  than  in 
this  example,  the  piezoelectric  wafers  are  used  to  monitor  deformation  in  the  thick-section 
sandwich  beam.  Relatively  small  magnitude  of  electrical  and  mechanical  inputs  are  experienced 
by  the  piezoelectric  wafers  since  the  thick-section  of  sandwich  beam  undergoes  relatively  small 
deformations.  The  piezoelectric  wafers  can  be  used  to  monitor  time-dependent  performances  in 
the  sandwich  composites.  The  smart  sandwich  beam  is  first  subjected  to  a  temperature  change 
from  room  temperature  to  80°C,  once  the  beam  is  under  isothermal  temperature  at  80°C,  a 
mechanical  load  is  prescribed  in  order  to  simulate  three  point  bending  deformation.  Figure  64 
(top  right)  shows  the  creep  responses  of  the  sandwich  composites  obtained  from  the  multi-scale 
model  and  experiments.  Good  correlations  are  observed.  The  corresponding  lateral  displacement 
and  axial  stress  fields  are  shown  in  Fig.  64  (bottom  left).  Finally,  the  recorded  electric  potentials 
in  the  piezoelectric  wafers  during  creep  deformation  are  also  reported. 
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Figure  64  A  multi-scale  model  for  predicting  time-temperature  dependent  response  of  smart  sandwich 
composites:  a  multi-scale  framework  (top  left);  overall  creep  response  in  smart  sandwich  beam  at  80“C 
(top  right);  overall  displacement  and  stress  fields  at  two  different  times  (bottom  left);  and  electric 
potential  reading  from  the  piezoelectric  sensor  during  the  creep  deformation  (bottom  right). 


b)  Active  fiber  composite 

The  nonlinear  time-dependent  electro-meehanieal  and  viseoelastic  eonstitutive  models,  whieh  are 
implemented  in  FE,  are  also  used  to  study  the  overall  eleetro-meehanieal  responses  of  AFCs, 
comprising  of  PZT  fibers,  epoxy  matrix,  and  electrode  fingers,  as  illustrated  in  Fig.  1.  The  AFC 
specimens  are  subjected  to  cyclic  electric  fields  with  different  amplitude  and  frequencies,  and  the 
overall  electro-mechanical  responses  are  monitored.  A  representative  unit-cell  model  of  AFCs  is 
considered,  as  shown  in  Fig.  65,  and  3D  continuum  elements  are  used  to  generate  the  unit-cell 
model  of  AFCs.  Electric  field  lines  are  obtained  through  the  positive  and  negative  electric 
potentials  applied  to  the  electrodes,  as  illustrated  in  Fig.  65.  The  nonlinear  time-dependent 
electro-mechanical  constitutive  model  (Section  4.2)  is  considered  for  the  PZT,  while  a 
viscoelastic  model  is  used  for  the  epoxy.  The  mechanical  responses  of  the  electrode  fingers  are 
modeled  as  linear  elastic,  while  a  linear  electro-static  model  is  considered  for  the  electrical 
behavior  of  the  electrode  fingers.  Figure  66  illustrates  the  hysteretic  responses  of  the  AFCs  from 
experiments  and  FE  analyses. 
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Figure  65  AFC  unit-cell 
Electromechanical  hysteresis  for  AFC  in  1Hz 
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Figure  66  Electro-mechanical  hysteretic  responses 


c)  Active  truss  systems 

The  finite  element  analyses  for  electro-active  materials  and  structures  are  also  presented  to  study 
shape  changes  in  smart  structures  comprising  of  3D  active  truss  systems.  Trusses  are  known  as 
load  bearing  components  in  many  structures.  Truss  systems  consist  of  relatively  slender 
members  connected  by  flexible  joints  with  rotational  and  translational  motions,  making  them 
appealing  for  compliant  structures.  Recent  advances  in  active  materials,  such  as  shape  memory 
and  electro-active  materials,  allow  for  generating  compliant  structures  that  can  sense  and  adapt  to 
various  stimuli  by  controllable  configurations  into  various  shapes.  There  have  been  several 
studies  on  understanding  the  geometry  of  the  truss  systems  for  certain  applications.  Sofia  et  al. 
(2009)  studied  morphing  hinged  truss  structures,  where  they  have  used  shape  memory  wires  in 
order  to  activate  their  structure.  Macareno  et  al.  (2008)  considered  a  linear  truss  made  of  3D 
tetrahedral  units  for  Variable  Geometry  Trusses  (VGTs).  Aviles  et  al.  (2000)  investigated  the 
position  problems  in  the  open-loop  variable  geometry  trusses.  An  optimization  scheme  is 
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designed  in  order  to  minimize  the  actuator's  displacement  and  consequently  its  energy.  This 
method  has  been  used  for  different  truss  architectures  in  specific  to  modular  tetrahedral  linear 
truss.  Reconfiguration  of  the  shape  of  structures  can  change  the  way  and  even  intensity  of  the 
load  that  the  structures  can  withhold.  Shape  changing  capability  can  improve  the  performance  of 
the  structures  and  is  expected  to  enhance  durability  of  the  structures  by  continuously  minimizing 
stresses  in  the  structures. 

This  study  focuses  on  theoretical  and  computational  methods  on  controlling  shape  changes  in 
compliant  truss  systems  comprising  of  electro-active  (piezoelectric)  materials.  An  elastic 
response  is  considered  for  the  mechanical  behavior  of  the  truss,  while  a  nonlinear  electro¬ 
mechanical  coupling  is  used  for  the  piezoelectric  components.  A  desired  shape  is  determined 
through  the  deformation  paths,  and  the  corresponding  strains  and  stresses  are  calculated.  The 
stresses  are  then  prescribed  by  applying  electric  field  though  the  piezoelectric  materials,  and  as  a 
result  the  truss  system  undergoes  the  desired  shape  changes.  The  procedure  is  summarized  as: 

1)  Define  equations  that  form  shape  changes  from  reference  configuration  (X,)  to  deformed 
configuration  (x,) 

1 

2)  The  deformation  gradient  and  strain  are  defined  as:  F  =  and  £  =  ~  l) 

3)  Unit  vector  direction  for  each  truss  element  is  determined  from  the  current  location 

4)  Determine  the  amount  and  direction  of  electric  field  that  should  be  prescribed  to  the  individual 
active  member  from  the  unit- vector  and  constitutive  relation  G(s,E)  =  0 

The  first  example  that  has  been  considered  in  this  study  is  a  linear  tetrahedral  truss.  The  desired 
shape  of  the  truss  is  considered  to  be  a  curve  that  is  clamped  at  the  origin  and  it  passes  the  points 
(L/4J./S),  (3Z/4,-Z/8),  (LJ./S)  with  L  being  the  total  length  of  the  linear  truss.  This  shape  is 
shown  in  Fig  61  (top).  This  truss  consists  of  70  tetrahedral  units,  72  nodes  and  213  truss 
members.  The  beam  theory  approximation  for  the  stresses  induced  in  the  beam  is  used  to  actuate 
shape  changing  in  the  truss,  as  shown  in  Fig.  67.  Another  example  is  done  on  generating  a 
sinusoidal  shape  change  by  electrically  activated  piezoelectric  components  in  the  truss  system,  as 
illustrated  in  Fig.  68.  It  is  seen  that  a  relatively  large  displacement  can  be  attained  while 
maintaining  a  relatively  small  strain  in  each  truss  component.  The  amount  of  voltage  needed  is 
based  on  the  electro-mechanical  properties  of  a  piezoelectric  tube  actuator. 


Figure  67  Beam  like  tmss:  Reference  configuration  (top)  and  deformed  configuration  (bottom) 
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Figure  68  Sinusoidal  shape  ehanges  produeed  by  aetuating  a  tetrahedral  truss  system 


d)  Pinching  of  a  slender  honeycomb  structure 

Another  example  of  inducing  lateral  deflections  (shape  changes)  in  a  slender  structure  is  by 
deforming  a  mid-section  of  a  slender  plate  or  beam  comprising  of  honeycomb-like  structures,  as 
illustrated  in  Fig.  69.  FE  analyses  of  smart  structures  are  used  to  study  the  lateral  deflection  of 
the  plate.  Consider  a  plate  with  actuators  placed  in  its  mid-section,  which  can  be  activated  in  the 
transverse  direction  of  the  plate.  The  honeycomb-like  microstructures  of  the  plate  results  in 
lightweight  and  compliant  systems.  By  deforming  the  mid-section  of  the  plate,  in  this  case 
through  the  piezoelectric  actuators,  the  plate  experiences  flapping  (bending).  The  amount  of 
lateral  deflection  at  the  free  end  that  can  be  gained  by  laterally  deforming  the  mid-section  is 
shown  in  Fig.  70.  It  is  seen  that  a  relatively  small  deformation  induced  on  the  mid-section  can 
generate  a  relatively  large  lateral  displacement.  Different  microstructural  designs  of  the  active 
slender  plate  can  also  be  considered  in  order  to  achieve  bending/flapping  of  the  plate.  One  of  the 
advantages  of  this  plate  design  is  that  it  is  relatively  easy  to  control  the  deformation  through 
electric  field  inputs  prescribed  at  one  location,  instead  of  applying  different  electric  field  inputs 
at  various  locations  in  order  to  achieve  desired  shapes,  like  in  the  active  truss  systems  or 
composite  beams  with  piezoelectric  patches. 
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Figure  69  A  slender  honeyeomb  strueture  undergoing  bending  due  to  pinebing 
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Figure  70  Lateral  deflection  obtained  from  pinebing  tbe  mid-point  of  tbe  beam 


4.5)  Concluding  Remarks 

Smart  composites  comprising  of  piezoelectric  materials  as  the  electro-active  component  are 
promising  for  multifunctional  systems  such  as  morphing  structures.  Lightweight  and  flexible 
active  structures  are  necessary  for  generating  morphing  systems  with  the  ability  to  respond  to 
various  external  stimuli.  Various  types  of  electro-active  lightweight  and  flexible  composites  have 
been  considered,  such  as  piezocomposites  comprising  of  piezoelectric  ceramics  inclusions  and 
ductile  matrix  (polymers  and  metals),  laminated  composites  with  integrated  piezoelectric 
materials,  arrangements  of  flexible  systems  with  active  components,  etc.  Several  experimental 
evidences  show  that  piezoelectric  ceramics  can  experience  nonlinear  and  time-dependent 
responses  when  subjected  to  electro-mechanical  responses.  High  electric  field  inputs  and 
compressive  stresses  can  lead  to  polarization  switching,  disrupting  the  electro-mechanical 
coupling  performance  of  the  multifunctional  systems.  Furthermore,  polymers  are  promising  for 
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generating  flexible  and  lightweight  systems;  however,  they  experienee  signifieant  viseoelastic 
responses.  Both  piezoeleetric  eeramies  and  polymers  dissipate  energy  when  subjeeted  to 
mechanical  and  non-mechanical  stimuli,  which  is  converted  into  heat  and  thus  increasing  the 
temperatures  of  the  systems.  It  is  also  known  that  the  electro-mechanical  performances  of 
piezoelectric  ceramics  and  polymers  change  with  temperatures.  The  differences  in  the  nonlinear 
time-dependent  thermo-electro-mechanical  responses  of  piezoelectric  ceramics  and  polymers, 
together  with  various  microstructural  morphologies  will  have  significant  effects  on  the  overall 
performance  of  multifunctional  morphing  systems. 

This  project  is  aimed  at  understanding  the  effect  of  coupled  thermal,  electrical,  and  mechanical 
responses,  including  loading  rate  (time)  effect,  of  the  constituents  (piezoelectric  ceramics  and 
polymers)  on  the  overall  multi-field  responses  of  electro-active  composites.  For  this  purpose,  a 
multi-scale  framework  that  consists  of  constitutive  models  for  constituents  incorporating  coupled 
deformation,  heat  conduction,  and  piezoelectric  analyses  with  stress,  temperature,  electric  field, 
and  time  dependent  responses;  nonlinear  micromechanical  models  for  various  composite 
reinforcements  including  functionally  graded  composites;  and  large  scale  multifunctional 
structural  analyses  has  been  formulated.  Limited  experiments  on  piezocomposites  have  been 
conducted  in  order  to  study  micro-scale  response  of  the  composite  and  provide  useful  insight  into 
the  effect  of  microstructural  morphologies  on  the  overall  (macroscopic)  multi-field  response  of 
composites. 

Phenomenological  models  within  continuum  mechanics  and  thermodynamics  approaches  have 
been  formulated  for  ferroelectric  ceramics  undergoing  various  levels  of  electric  field  inputs  and 
mechanical  stresses,  and  for  polymers  under  various  histories  of  mechanical  loadings.  For  the 
ferroelectric  ceramics,  the  responses  from  the  constitutive  models  have  been  successfully 
compared  with  available  experimental  data,  incorporating  time  (rate),  temperature,  and  electric 
field  dependent  behaviors.  Fully  coupled  thermo-viscoelastic  models  have  been  formulated,  and 
shown  capable  in  capturing  the  overall  heat  conduction  and  deformation  in  polymers  and 
polymeric  composites.  Furthermore,  efficient  numerical  integration  algorithms  have  been 
developed  in  order  to  simulate  the  overall  performance  of  ferroelectric  ceramics  and  polymers 
under  various  histories  of  external  stimuli.  The  numerical  integration  algorithms  are  easily 
implemented  within  finite  element  and  finite  difference  numerical  methods,  and  integrated  to 
various  micromechanics  models. 

Several  micromechanics  models,  i.e.,  simplified  unit-cell  methods  and  micromechanics  models 
with  microstructural  details,  have  been  considered  in  order  to  determine  the  overall  nonlinear 
time-dependent  behaviors  of  composites  undergoing  couple  heat  conduction  and  electro¬ 
mechanical  deformations.  The  micromechanical  models  incorporate  different  nonlinear  time-  and 
field  dependent  constitutive  models  for  the  inclusions  and  matrix.  Efficient  numerical  algorithms 
have  been  formulated  for  different  micromechanical  models,  and  implemented  in  finite  element 
analyses.  The  effective  (macroscopic)  responses  determined  from  the  unit-cell  models  have  been 
compared  with  the  overall  responses  from  available  experimental  data  and  micromechanics 
models  with  microstructural  details.  The  unit-cell  models  are  shown  capable  in  capturing  the 
macroscopic  responses  of  the  composites  under  various  histories  of  loadings,  even  though 
relatively  crude  approximations  are  considered  in  formulating  the  unit-cell  models.  The 
micromechanics  studies  also  reveal  the  importance  of  thermal  stress  effects,  which  occurs  due  to 
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temperature  changes  during  either  processing  of  the  composites  or  in-service  conditions,  on  the 
overall  electro-mechanical  responses  of  composites.  It  is  noted  that  relatively  high  thermal 
stresses  in  the  constituents  of  the  composites  often  lead  to  significant  nonlinear  responses, 
including  possibilities  of  damage/degradation,  and  accelerate  the  time-dependent  deformations, 
and  thus  affecting  the  overall  responses  of  the  composites.  The  unit-cell  models,  however,  are 
limited  in  capturing  the  discontinuities  of  field  variables  at  the  interfaces  between  the  inclusions 
and  matrix,  and  localized  field  variables,  such  as  stress  concentrations.  Discontinuities  and 
localized  field  variables  often  cause  debonding  and  damage  of  the  composites.  One  of  the  major 
advantages  in  using  the  simplified  unit-cell  models  over  the  micromechanical  models  with 
microstructural  details  is  in  cutting  the  computational  cost,  while  reasonably  predicting  the 
overall  nonlinear  and  time-dependent  responses  of  the  composites.  Thus,  the  unit-cell  models  are 
promising  for  analyzing  large-scale  structures  under  various  histories  of  mechanical  and  non¬ 
mechanical  stimuli,  and  at  the  same  time  incorporating  the  nonlinear  and  time-dependent 
responses  of  the  constituents. 

A  multi-scale  framework  has  been  formulated  by  integrating  nonlinear  time-dependent 
constitutive  models,  micromechanics  models,  and  structural  analyses  via  finite  element  and  finite 
difference  methods.  The  multi-scale  framework  is  used  to  study  the  overall  performance  of 
active  structures  or  structural  components  undergoing  various  histories  of  external  stimuli,  and 
perform  active  controls  using  the  active  materials,  i.e.,  piezoelectric  ceramics  and  piezoelectric 
composites,  in  order  to  achieve  desired  performance  of  the  structures.  Several  examples  have 
been  presented  on  the  ability  of  the  multi-scale  framework  in  simulating  time-dependent  multi- 
field  responses  of  active  structures  or  structural  components.  It  is  noted  that  for  precisely 
controlling  the  deformed  shapes  of  structural  components  using  electro-active  composites,  it  is 
necessary  to  take  into  account  the  different  time-dependent  behaviors  of  the  constituents  and 
microstructural  morphologies  in  the  active  composites  since  they  significantly  affect  the  overall 
performance  of  the  smart  structures. 
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